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1 .0  Before  You  Start 


BIOSCREEN  is  a  software  tool  to  be  used  to  model  groundwater  plume  migration,  taking  into 
account  the  effects  of  natural  attenuation  and  biodegradation  of  mobile  contaminants.  The 
BIOSCREEN  system  consists  of  the  following  files: 


README.WRI  (this  file) 
BIOSCRN  XLS 
BIOSCRN  HLP 
BIOSCRN  PDF 
EXAMPLES  PDF 


Background  information  and  installation  instructions. 

The  BIOSCREEN  program  file,  in  Microsoft  Excel  5  0  format 
The  online  help  file,  in  Windows  Help  format. 

The  full  documentation,  in  Adobe  PDF  format 

Case  Study  information  (Appendix  A. 6)  in  Adobe  PDF  format 
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Run  SETUP  EXE  from  the  floppy  drive,  either  by  selecting  Run  from  the  File  menu  in  Program 
Manager  or  by  double-clicking  on  the  file  SETUP  EXE  in  File  Managr  r  (or  Windows  95  Explorer) 
The  installation  process  creates  the  C  \t5k  rSCRN  ■'ubdi rectory  on  your  hard  drive,  unless  you 
install  It  elsev/here.  and  copies  BIOSCRN  XLS,  BKrSCRN.HLP,  and  this  README.WRI  file  iniu 
ttie  new  diiectony 


To  run  RIOSCREEN  after  installation,  start  Microsoft  fExcol  and  open  the  BIOSCRN  XI  S  file  Fill 
in  your  field  or  hypothetical  data  into  the  blanks  provided  and  select  either  the  Run  Centerline 
(centerline  model)  or  Run  Array  (three-dimensional  plume  migration  model)  buttons  Loading  from 


x>iahv ' 

App»  0,1  'H'  r 
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within  Excel  is  preferred  to  double-clicking  the  filename  in  file  manager  (see  Section  3.0  for  an 
explanation  of  the  evils  you  v.'ill  avoid).  Alternatively,  you  can  add  your  new  BIOSC  REEN 
directory  to  your  path. 

1 .2  Minimum  System  Requirements 

Any  PC-Compatible  computer  system  capable  of  running  Microsoft  Excel  version  5  0  for  Windows 
will  be  able  to  run  BIOSCREEN  Because  of  the  intensive  floating  point  calculations  used  in  the 
model,  however,  GSI  recommends  that  the  following  minimum  requirements  be  observed  to 
ensure  a  minimum  standard  of  performance: 

An  Iruel  486  CPU,  operating  at  60  MHz  or  faster 

8MB  Random  Access  Memory  (16  MB  or  more  will  improve  performance) 

1MB  free  hard  drive  space  (for  BIOSCREEN  only  -  Excel  will  require  substantially  more) 

1.3  Software  Requirements 

BIOSCREEN  is  designed  to  be  run  under  Microsoft  Excel  for  Windows.  Your  system  must  be 
running  Windows  3.1  or  Windows  for  Workgroups  3  11  and  have  Excel  for  Windows  version  5  0  or 
later  property  installed. 

BIOSCREEN  has  been  tested  on  Excel  version  5  0  for  Windows  on  systems  running  Windows 
3.1,  Windows  for  Workgroups  3.1 1,  Windows  95,  and  Windows  NT  Version  3.31  and  Excel  7  0  on 
Windows  95  and  Windows  NT  3.51.  Although  no  substantial  problems  were  seen  with  any  of 
these  configurations,  preliminary  reports  suggest  that  users  of  Excel  versions  prior  to  the  'c' 
maintenance  release  (versions  5  0  and  5  Oa)  may  experience  difficulties  under  some 
circumstances.  We  are  researching  this  issue. 


2.0  Installation 

1 .  Run  SETUP, EXE  from  the  Windows  Program  Manager  by  selecting  Run  from  the  File 
Menu  and  typing  A. SETUP,  or  run  it  from  File  Manager  by  double-clicking  its  icon 
This  will  create  the  directory  BIOSCRN  and  copy  the  essential  files  (listed  above)  to 
that  directory 

2  Start  Microsoft  Excel  and  open  the  BIOSCRN  XL8  ‘ile.  Fill  in  your  field  or  hypothetical 
data  into  the  blanks  provided  and  select  either  the  P,i.:n  Centerline  (centerline  2D 
model)  or  Run  Array  (three-dimensional  plume  nigi;;tion  model)  buttons.  The 
unshaded  (white)  cells  are  intended  to  receive  your  data,  while  light  grey  cells  may 
contain  either  your  data  or  a  default  formula  i,r«-iii sorted  into  the  cells  by  clicking  the 
"Restore  Default  Formulas"  button  on  the  main  input  screen)  Dark  grey  cells  contain 
only  formulas  used  by  the  BIOSCREEN  system,  and  should  not  be  changed  by  the 
user. 

3  You  can  load  and  view  the  online  Help  by  double-clicking  its  fit  mame  in  the  File 
Manager  pane  or  by  selecting  the  help  button  fiom  the  BIOSCREEN  inteiface 


3  0  Troubleshooting 

When  I  tty  to  get  online  help,  the  Windows  help  program  starts,  but  a  dialog  box  tells  me  it 
cannot  open  the  help  file.  Close  BIOSCREEN,  leaving  Excel  running.  From  the  File  menu, 
select  Open  Browse  the  directory  structure  until  you  find  BIOSCRN  XLS,  and  open  it.  Nov/,  the 
Help  button  will  open  the  proper  help  file.  See  online  help  for  additional  (more  permaneni) 
met'  ’ds  of  solving  this  problem. 

The  Input  screen  appears,  but  it  Is  very  small  or  its  edges  extend  off  the  screen.  The 

system  was  designed  to  operate  at  a  resolution  of  640x480  pixels,  or  standard  VGA  resolution  If 
you  are  using  a  higher  resolution,  you  may  either  change  your  video  driver  m  Windows  Setup  to 
standard  VGA  resolution  or  modify  the  ZoomFactor  ‘o  accci  it  foi  your  configuration  See  online 
tielp  for  Display  configuration  information  and  instructions  for  changing  the  ZoomFac.or 
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Additional  Information 
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4.1  Printing  the  BIOSCREEN  Data  Screens 

If  you  want  to  print  the  data  in  your  input  screen,  the  centerline  model  screen,  or  the  3-D  graph 
screen,  a  named  range  for  each  has  been  pre-configured.  By  pressing  r5  and  selecting  the 
Print__Area  range,  the  screen  you  are  printing  will  be  highlighted.  The  same  range  name  exists  for 
each  screen.  Printing  any  other  part  of  the  worksheets  will  require  resetting  the  print  area 
Consult  the  Excel  documentation  for  instructions  on  how  to  accomplish  this. 

4.2  Electronic  Manual  Availability 

The  full  documentation  set  for  BIOSCREEN  is  also  available  rn  Adobe  Portable  Document  Format 
(PDF).  Two  separate  files  are  available,  the  manual  and  a  set  of  two  sample  problems  which 
illustrate  a  realistic  application  of  BIOSCREEN.  Consult  the  on-line  help  file  (BIOSCRN  HLP)  for 
sources  for  this  electronic  documentation  Those  files  duplicate  the  contents  of  the  printed 
documentation. 
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INTRODUCTION 


BIOSCREEN  is  an  easy-tcv-use  screening  model  which  simulates  remediation  through  natural 
attenuation  (RNA)  of  dissolved  hydrocarbons  at  petroleum  fuel  release  sites.  The  softA-are, 
programmeii  in  the  Microsoft*  Excel  spreadsheet  environment  a^id  based  on  the  Domenico 
analytical  solute  transport  model,  has  the  ability  to  simulate  advection,  dispersion, 
adsorption,  and  aerobic  decay  as  well  as  anaerobic  reactions  that  have  been  shown  to  be  the 
dominant  biodegi  'dation  processes  at  many  petroleum  release  sites.  BIOSCREEN  includes 
three  different  model  types: 

31  Solute  transport  without  decay. 

21  Solute  transport  with  biodegradation  modeled  as  a  first-order  decay  process  (simple,  lumped-parameter 
approach), 

3)  Solute  transport  with  biode'eradatioii  modeled  as  an  "instantaneous'  biodegradation  reaction  (approach 
used  by  BKjPlUML  modets) 

The  model  is  designed  to  simulate  biodegradation  by  both  aerobic  and  anaerobic  reactions.  1 1 
was  developed  for  the  Air  Force  Center  for  Environmental  Excellence  (AFCEE)  Technology 
Transfer  Division  at  Brooks  Air  Force  Base  by  Groundwater  Services,  Inc.,  Houston,  Texas. 

INTENDED  USES  FOR  BIOSCREEN 

BIOSCREEN  attempts  to  answer  two  fundamental  questions  regarding  RNA: 

1 .  How  far  will  the  dissolved  contaminant  plume  extend  if  no  engineered  controls 
or  further  source  rone  reduction  measures  are  implemented? 

BIOSCREEN  use.s  an  analytical  solute  transport  model  with  two  options  for  simulating 
in-situ  biodegradation:  first-order  decay  and  instantaneous  reaction.  The  model  will 
predict  the  maximum  extent  of  plume  migration,  which  may  then  be  compared  to  the 
dist.mce  to  potential  points  of  exposure  (e.g.,  drinking  water  wells,  groundwater 
di.scharge  areas,  or  property  boundaries).  Analytical  groundwater  transport  models 
have  .seen  wide  application  for  this  purpo.st*  (e  g.,  ASTM  19*^5),  and  experience  has 
shown  such  models  can  produce  reliable  results  when  site  conditions  in  the  plume  area 
are  relatively  uniform. 

2.  How  long  will  the  plume  persist  until  natural  attenuation  processes  cause  it  to 
dissipate? 

BIOSCREEN  uses  a  simple  ma.s.s  balance  approach  bn.sed  on  the  mass  of  dissolvable 
hydrocarbons  in  the  source  zone  and  the  rate  of  hydrocarbons  le  iving  the  source  zone  to 
estimate  the  source  zone  concentration  vs.  time.  Because  an  exponential  decay  in  sounre 
zone  concentration  is  assumed,  the  predicted  plume  lifetimes  can  be  large,  usually 
rangmg  from  5  to  500  years.  Note  this  is  an  unverified  relationship  as  there  are  few 
data  showing  source  concentrations  vs.  long  time  periods,  and  tlu-  lesulis  should  be 
considered  order  of-magnitude  estimates  of  the  time  n  mired  to  dissipate  the  plume. 

BIOSCREEN  is  intended  to  be  used  m  two  wavs: 

I .  As  a  screening  model  to  determine  if  RNA  is  feasible  at  a  site. 

In  thi.s  case,  BiOSCREEN  is  used  early  in  the  remedial  investigation  to  determine  if  an 
RNA  field  program  should  be  implemented  to  quantify  the  natural  attenuation 


ocairring  at  a  site.  Some  data,  such  as  electron  acceptor  concentrations,  may  not  be 
available,  so  typical  values  are  used.  In  addition,  the  model  can  be  used  help 
develop  long-term  monitoring  plans  for  RNA  projects. 

2.  As  the  primary  RNA  groundwater  model  at  smaller  sites. 

The  Air  Force  Intrinsic  Remediation  Protocol  (Wiedemeier,  Wilson,  et  a!.,  1995) 
describes  how  groundwaier  models  may  be  used  to  help  verify  that  natural  attenuation 
is  occurring  and  to  help  predict  how  far  plumes  might  extend  under  an  RNA  scenario. 
At  large,  high-effort  sites  such  as  Superfund  and  RCRA  sites,  a  more  sophisticated 
model  such  as  BIOPLUME  is  probably  more  approprir  j.  At  less  complicated,  lower- 
effort  sites  such  as  service  stations,  BIOSCREEN  may  be  sufficient  to  complete  the 
RNA  study.  (Note;  "Intrinsic  remediation"  is  a  risk-based  strategy  that  relies  cn 
RNA). 

BIOSCREEN  ha.s  the  followmg  limitations: 

1.  As  an  analytical  model,  BIOSCREEN  a.ssumes  simple  groundwater  flow 
conditions. 

The  model  should  not  be  applied  where  pumping  systems  create  a  complicated  flow 
field.  In  addition,  the  model  should  not  be  applied  where  vertical  flow  gradients 
aftect  contaminant  transport. 

2.  As  an  screening  tool,  BIOSCREEN  only  approximates  more  complicated 
processes  that  occur  in  the  field. 

Tlie  mtxlel  should  not  be  applied  where  extremely  detailed,  accurate  results  that 
closely  match  site  conditions  arc  requiied.  More  comprehensive  numerical  moilels 
should  be  applied  in  these  cases. 


FUNDAMENTALS  OF  NATURAL  ATTENUATION 

Biodegradation  Modeling 

Naturally  occurring  biological  pnxes.ses  can  significantly  enhance  the  rate  of  organic  mass 
removal  from  contammated  aquifers.  Biodegrad.ition  research  performed  by  Rice  Universit), 
government  agencies,  and  other  research  groups  has  dentified  several  n' iin  themes  that  are 
crucial  for  future  studies  of  natural  attenuation: 

1.  The  relativt'  impcrtancc  of  ynnindwater  transport  vs.  microbial  kinetics  is  a  key 
consideration  for  developing  workable  biodegradation  expressions  in  models.  Results 
from  the  United  Creosote  site  (Texas)  and  the  Traverse  City  Tiiel  Spill  site  (Michigan) 
indicate  that  biodegradation  is  better  represented  as  a  macro-scale  wastewater 
treatment-type  vroress  Hian  as  a  micro-scale  study  of  microbial  reactions. 

2.  rill  distribution  and  availability  of  electron  acceptors  control  the  rate  of  in-situ 
biodegradation  for  most  petroleum  release  site  plumes.  i)t!icr  factors  (e.g.,  papula tioii 
of  microbes,  pH,  temperature,  etc.)  rarely  limit  the  amount  of  biodegradati'm  occurring 
at  these  sites. 

Borden  et  al.  (NS6)  developied  the  BIOPLUME  model,  which  simulate.s  aerobic  biodegradation 
as  an  "instantaneous"  microbial  reaction  that  i.s  limited  by  the  .imouiat  of  electron  acieptor. 
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oxygen,  that  is  available.  In  other  words,  the  microbial  reaction  is  assumed  to  occur  at  a  much 
faster  rate  than  the  time  required  for  the  aquifer  to  replenish  tlie  amount  of  oxygen  in  the 
plume.  Although  the  time  required  for  the  biomass  to  aerobically  degrade  the  dissolved 
hydrocarbons  is  on  the  order  of  days,  the  overall  time  to  flush  a  plume  with  fresh  groundwater 
i.s  on  the  order  of  years  or  tens  of  years.  Borden  et  al.  (198b)  incorporated  a  simplifying 
a.ssumption  that  the  microbial  kinetics  are  instantaneous  into  the  USGS  two-dimensional  solute 
transport  model  (Konikow  and  Bredehoeft,  1978)  using  a  simple  superposition  algorithm.  The 
resulting  model,  BIOPLUME,  was  able  to  simulate  solute  transport  and  fate  under  the  effects  of 
instantaneous,  oxygen-limited  in-situ  biodegradation. 

Rifai  and  Bedient  (1990)  extended  this  approach  and  developed  the  BIOPLUME  II  model, 
which  simulates  the  transport  of  two  plumes:  an  oxygen  plurrie  and  a  contammant  plume.  The 
two  plumes  are  allowed  to  react,  and  the  ratio  of  oxygen  to  contaminant  consumed  by  the 
reaction  is  determined  from  an  appropriate  stoichiometric  model.  The  BIOPLUME  II  mixfel  is 
documented  with  a  detailed  user's  manual  (Rifai  et  al.,  1987)  and  is  currently  being  used  by  EPA 
regional  offices,  U.S.  Air  Force  facilities,  and  bv  consulting  firms.  Borden  et  al.  (1986)  applied 
the  BIOPLUME  concepts  to  the  Conroe  Superfund  site;  Rifai  et  nl.  (1988)  and  Rifai  et  at.  (1991) 
applied  the  BIOPLUME  II  model  to  a  jet  fuel  spill  at  a  Coast  Guard  facility  m  Michigan. 
Many  other  studies  using  the  BIOPLUME  II  model  have  been  presented  in  recent  literature. 

The  BIOPLUME  II  mtxlcl  has  increased  the  understanding  of  biodegradation  and  natural 
attenuation  by  simulating  the  effects  of  adsorption,  dispersion,  and  aerobic  biodegradation 
processes  m  one  model.  It  incorporates  a  simplified  mechanism  (first-order  decay)  for  handling 
other  degradation  processes,  but  does  not  address  specific  anaerobic  decay  reactions.  Early 
conceptual  models  of  natural  attenuation  were  based  ai  the  assun'iption  that  the  anaerobic 
degradation  pathways  were  too  slow  to  have  any  meanmgful  effect  on  the  overall  natural 
attenuation  rate  at  most  sites.  Accordingly,  most  field  programs  focused  only  on  the  distribution 
of  oxygen  and  contaminants,  and  did  not  measure  the  indicators  of  anaerobic  activity  sudi  as 
depletion  of  anaerobic  electron  acceptors  or  accumulation  of  anaerobic  meMbolic  by-products. 

The  Air  Force  Natural  Attenuatiun  Initiative 

Gver  the  past  several  years,  the  high  cost  and  poor  performance  of  many  pump-and-treat 
remediation  systems  have  led  many  researchers  to  consider  RNA  as  an  alternative  technology 
tor  groundwater  remediation.  A  detailed  understanding  of  natural  attenuation  processes  is 
needed  to  support  the  development  of  this  remediation  approach.  Researchers  associated  with 
the  U.S.  EPA's  R.S.  Kerr  Environmental  Research  Laboratory  (now  the  Subsurface  Protection 
and  Remediation  Division  of  the  National  Risk  Management  Laboratory)  have  suggested  that 
anaerobic  pathways  could  be  a  significant,  or  even  the  dominant,  degradation  mechanism  at 
many  petroleum  fuel  sites  (Wilson,  1994).  The  natural  attenuation  initiative,  developed  by  the 
AFCEE  Technology  Transfer  Division,  was  designed  to  investigate  how  natural  attenuation 
processes  affect  the  migration  of  plumes  at  petroleum  release  sites.  Under  the  guidance  of  l.t. 
Col.  Ross  Miller,  a  tliree-pronged  technology  development  effort  was  launched  in  1993  which 
will  ultimately  consist  of  the  following  elements: 

1)  Field  data  collected  at  over  30  sites  around  the  country  (Wiedemeier,  Miller,  et  al., 
1 99.5  >  analyzinf^  aerobic  and  anaerobic  processes. 

2)  A  Technical  Protocol,  outlining  the  approach,  data  collection  techniques,  and  data 
aualusis  methods  required  for  conducting  an  Air  Force  RNA  Utiidy  ( W icd 'tneicr , 
Wilson,  et  al.,  19951 

■  JJ  lioo  RNA  modeling  tools,  the  BIOPLUMF  III  model  beinq  developed  by  Dr.  lianadi 
Rifai  al  Rice  Universitq  (Ril.u  et  al.,  199.5),  and  the  BIDSCRF.F.N  model  developed  by 
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Groundwater  Services,  Inc.  (BIOPLUME  III,  a  more  sophisticated  biodegradation 
model  than  BIOSCREEN,  employs  particle  tracking  of  both  hydrocarbon  and  alternate 
electron  acceptors  using  a  numerical  solver.  The  model  employs  sequential  degradation 
of  the  biodegradation  reactions  based  on  zero  order,  first  order,  instantaneous,  or  Monod 
kinetics). 

ReUtive  Importance  of  Different  Electron  Acceptors 

The  Intrinsic  Remediation  Technical  Protocol  and  modeling  tools  focus  on  evaluating  both 
aerobic  (in  the  presence  of  oxygen)  and  anaerobic  (without  oxygen)  degradation  processes.  In 
the  presence  of  organic  substrate  and  dissolved  oxygen,  microorganisms  capable  of  aerobic 
metabolism  will  predominate  over  anaerobic  forms.  However,  dissolved  oxygen  is  rapidly 
consumed  in  the  interior  of  contaminant  plumes,  converting  these  areas  into  anoxic  (low-oxygen) 
zones.  Under  these  conditions,  anaerobic  bacteria  begin  to  utilize  other  electron  acceptors  to 
metabolize  dissolved  hydrocarbons.  The  principal  factors  influencing  the  utilization  of  the 
various  electron  acceptors  by  fuel-hydrocarbon-degrading  bacteria  include-  1)  the  relative 
biochemical  energy  provided  bv  the  reaction,  2)  the  availabilitv  of  individual  or  specific 
electron  acceptors  at  a  particular  site,  tind  3)  the  kinetics  (rate)  of  the  microbial  reaction 
associated  with  the  different  electron  acceptors. 

Preferred  Reaction.':  by  Energy  Potential 

Biologically  mediated  degradation  reactions  are  reduction/oxidation  (redox)  reactions, 
involvmg  the  transfer  of  electrons  from  the  orgaruc  contaminant  compound  to  an  electron 
acceptor.  Oxygen  is  the  electron  acceptor  for  aerobic  metabolism,  whereas  nitrate,  ferric  iron, 
sulfate,  and  carbon  dioxide  can  serve  as  electron  acceptors  for  .alternative  anaerobic  pathways. 
This  transfer  of  electrons  releases  energy  which  is  utilized  for  microbial  cell  maintenance  and 
growth.  The  biochemical  energy  associated  with  alternative  degradation  pathways  can  lx 
represented  by  the  redox  potential  of  the  alternative  electron  acceptors:  the  more  positive  the 
redox  potential,  tile  more  energetically  favorable  the  re.iction.  With  everything  else  being 
equal,  organisms  with  more  efficient  modes  of  metabolism  grow  faster  and  therefore  dominate 
over  less  efficient  forms. 


Electron 

Acceptor 

Type  of 
Reaction 

Metabolic 

By-Product 

Redox  Potential 
(pH  =  7,  In  mvolta)* 

Reaction 

Preference 

Oxygen 

Aerobic 

COn 

+  820 

Most  Preferred 

Nitrate 

Anaerobic 

N2,  UOn 

+  740 

11 

Ferric  Iron 
(solid) 

Anaorobic 

Ft’rrims  Iron 
(dissolved) 

-50 

U 

Sulfate 

Anaenibic 

lUS 

-  220 

U 

Carbon  Dioxide 

Anaerobic 

Methane 

-  240 

Least  Preferred 

*  from  W,edemeier,  Wilson,  et  al.,  1995. 


Based  .solely  cn  thermodynamic  considerations,  the  most  energetically  preferred  reaction 
should  prrx'eed  in  the  plume  until  all  of  the  required  electron  acceptor  is  depleted.  At  that 
pomt,  tlie  next  most-preferred  reaction  should  begin  and  continue  until  that  electron  acceptor  is 
consumed,  leading  to  a  pattern  where  preferred  eleclroii  acceptors  are  consumed  one  at  a  time,  in 
sixjuence.  Based  in  thi.s  principle,  one  would  expect  to  observe  monitoring  well  data  with  no 
detect "  results  tor  the  more  energetic  electron  acceptors,  such  as  oxygen  and  nitrate,  in  locations 
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where  evidence  of  less  energetic  reactions  is  observed  (e.g.  monitoring  well  data  indicating  the 
presence  of  ferrous  iron). 

In  practice,  however,  it  is  unusual  to  collect  samples  from  monitoring  wells  that  are  completely 
depleted  in  one  or  more  electron  acceptors.  Two  processes  are  probably  responsible  for  this 
observation; 


L  Alternative  biochemical  mechanisms  exhibiting  very  similar  energy  potentials  (such 
as  aerobic  oxidation  and  nitrate  reduction)  may  occur  concurrently  when  the  preferred 
electron  acceptor  is  reduced  in  concentration,  rather  than  fully  depleted.  Facultative 
aerobes  (bacteria  able  to  utilize  electron  acceptors  in  both  aerobic  and  anaerobic 
environments),  for  example,  can  shift  from  aerobic  metabolism  to  nitrate  •■eduction 
when  oxygen  is  still  present  but  at  low  concentrations  (i.e.  1  mg/L  oxygen:  Snoeyink  and 
Jenkins,  1980).  Similarly,  noting  the  nearly  equivalent  redox  potentials  for  sulfate  and 
carbon  dioxide  (-220  millivolts  and  -240  millivolts,  respectively)  one  might  expect  that 
sulfate  reduction  and  methanogenic  reactions  may  also  occur  together. 


2.  Standard  monitoring  wells,  with  5-  to  10-  foot  screened  intervals,  will  mix  waters  from 
different  vertical  zones.  If  different  biodegradation  reactions  are  occurring  at  different 
depths,  then  one  would  expect  to  find  geochemical  evidence  of  alternative  degradation 
mechanisms  occurring  in  the  same  well.  If  the  dissolved  hydrocarbon  plume  is  thinner 
than  the  screened  interval  of  a  monitoring  well,  then  the  geochemical  evidence  of 
electron  acceptor  depiction  or  metabolite  accumulation  will  be  diluted  by  mixing  with 
clean  water  from  zones  where  no  degradation  is  occurring. 


Therefore,  most  natural  attenuation  programs  yield  data  that  indicate  a  general  pattern  of 
electron  acceptor  depletion,  but  not  complete  depletion,  and  an  overlapping  of  electron 
acceptor/metabolite  isopleths  into  zones  not  predicted  by  thermodynamic  principles.  For 
example,  a  zone  of  methane  accumulahon  may  be  larger  than  the  apparent  .moxic  zone. 
Nevertheless,  these  general  patterns  of  geochemical  changes  within  the  plume  area  provide 
strong  evidence  that  multiple  mechanisms  of  biodegradation  are  occurring  at  many  sites.  The 
BIOSCREEN  software  attempts  to  account  for  the  majority  of  these  biodegradation 
mechanisms. 


Distribution  of  Electron  Acceptors  at  Sites 

The  utilization  of  electron  acceptors  is  generally  based  on  the  energy  of  the  reaction  iind  the 
availability  of  the  electron  acceptor  at  tlie  site.  While  the  energy  of  each  reaction  i.s  based  in 
thermodyi  imics,  the  distribution  of  electron  acceptors  is  dependent  cn  site-specific 
hydrogeochemical  processes  and  can  vary  significantly  among  sites.  For  example,  a  study  of 
several  sites  yielded  the  following  summary  of  available  electron  acceptors  and  metabolic  by¬ 
products; 
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Measured  Background  Electron  Acceptor/By-Product  Concentration  (mg/t) 

Base  Facility 

Background 

Oxygen 

Background 

Nitrate 

Maximum 
Ferrous  Iron 

Background 

Sulfate 

Maximum 

Methane 

POL  Site, 

Hill  AFB,  Utah* 

6.0 

36.2 

55.6 

96.6 

2.0 

Hangar  10  Site, 

Elmendorf  AFB,  Alaska* 

0.8 

64.7 

8.9 

25.1 

9.0 

Site  ST-41. 

Elmendorf  AFB, Alaska* 

12. '7 

60.3 

40.5 

57.0 

\3 

Sue  ST-29. 

Patrick  AHb.  Florida* 

3.3 

0 

2.0 

0 

13.6 

Bldg.  735. 

Grissom  AFB.  Indiana 

9.1 

1.0 

2.2 

59.8 

1.0 

SW  MU  66  Site, 

Keesler  AFB,  MS 

1.7 

0.7 

,36.2 

22.4 

7.4 

POL  B  Site. 

7'yndall  AFB,  Florida 

1.4 

0.1 

13 

5.9 

4.6 

*Daa  collected  by  Pardons  Enginecnng  Science.  Inc.:  all  other  data  collected  by  Groundwater  Services.  Inc. 


At  tho  I’atrick  AFB  site,  nitrate  and  sulfate  are  not  important  electron  acceptors  while  the 
oxygen  and  the  methanogenii  reactions  dominate  (Wicdemeici",  Swanson,  etal.,  1995).  At  Hill 
AFB  cuid  CrLs.soni  AFB,  the  sullate  reactions  are  extremely  important  because  of  the  large 
amount  of  available  sulfate  for  reduction.  Note  that  different  sites  in  close  proximity  can  have 
quite  different  electron  acceptor  concentrations,  as  shown  by  the  two  sites  at  F.lmendorf  AFB. 
For  data  on  more  sites,  see  Table  1. 


Kinetics  of  Aerobic  and  Anaerobic  Reactions 

As  described  above,  aerobic  biodegradation  can  be  .simulated  as  an  'instantaneous"  reaciiuu 
that  is  limited  by  the  amount  of  electron  acceptor  (oxygen)  that  is  available.  Tne  microbial 
reaction  is  assumed  to  ixxur  at  a  much  faster  rate  than  the  time  required  for  the  aquifer  to 
replenish  the  amount  of  oxygen  in  the  plume  (Wilson  etal. ,  1985).  Although  the  time  required 
for  the  biomass  to  aerobically  degrade  the  dissolved  hydrocarbons  is  cn  the  order  of  days,  the 
overall  time  to  flush  a  plume  with  fresh  groundwater  is  on  the  order  of  years  or  tens  of  years. 

For  example,  micnxrosm  data  presented  by  Davis  ct  al.  (1994)  show  that  microbes  in  an 
environment  with  an  excess  of  electron  acceptors  c.m  degrade  high  concentration.s  of  dissolved 
heo/ene  very  rapidly  In  file  pre.senee  of  suqilus  oxygen,  aerobic'  bacteria  eon  degrade  ~1  nig/I, 
dissuh  ed  beiv.ene  m  about  h  davs,  whicli  can  be  considered  relatively  fast  (referred  to  as 
"instantaneous")  eompareti  to  the  years  required  for  flowing  groundwater  to  replenish  tin 
plume  area  with  oxygen. 


6 


•  •  • 


•  •  •  •  • 


L  tilizjfion  ficJars  nf  the  electron  icceptors/by-pnxJucts  are  (olloHS  I  mg  of  flev'triin  *ccrptt>r  Of  byproduct/mg  BTBX)  Dissolved  0»yger«  3  H.  Nitrate  4  9.  Iron  218,$u]fA'e  4  7,  Methane.  0  78 
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Recent  results  from  the  AFCEE  Natural  Attenuation  Initiative  indicate  that  tfie  anaerobic 
reactions,  which  were  originally  thought  to  be  too  slow  to  be  of  significance  in  groundwater,  can 
also  be  simulated  as  mstantsmeous  reactions  (Newell  et  al.,  1995).  For  example,  Davis  et  al. 
(1994)  also  ran  microcosm  studies  with  sulfate  reducers  and  methanogens  that  indicated  that 
benzene  could  be  degraded  in  a  peried  of  a  few  weeks  (after  acclimation).  When  compared  to 
the  time  required  to  replenish  electron  acceptors  in  a  plume,  it  appears  appropriate  to  simulate 
anaerobic  biodegradation  of  dissolved  hydrocarbons  with  an  instantaneous  i  action,  just  as  for 
aerobic  biodegradation  processes. 

This  conclusion  is  supported  by  observmg  the  pattern  of  anaerobic  electron  acceptors  and 
metabolic  by-products  along  the  plume  at  RNA  research  sites: 


If  microbiati  kinetics  were  limiting  the 
rate  of  biodegradation: 

•  Anaerobic  electron  acceptors  (nitrate  and 
sulfate)  would  be  constantly  decreasing  in 
concentration  as  one  moved  downgradient 
from  the  source  zone,  and 


If  microbial  kinetics  were  relatively  fast 
(instantaneous): 

•  Anaerobic  electron  acceptors  (nitrate  and 
sulfate)  would  be  mostly  or  totally 
consul!  led  in  the  source  zone,  and 


•  Anaerobic  by-products  (ferrous  iron  and 
methane)  would  be  constantly  increasing 
in  concentration  as  one  moved 
downgradient  from  the  source  zone. 


X 


•  Anaerobic  by-products  (ferrous  iron  and 
methane)  would  be  found  in  the  highest 
concentrations  in  the  source  zone. 


X - ► 


Tlie  second  pattern  is  observed  at  RNA  demon.stration  sites  (see  Figure  1),  supporting  the 
hypothesis  that  anaerobic  reactions  can  be  considered  to  be  relatively  instantaneous  at  most  or 
dimo.st  all  petroleum  release  sites.  From  a  theoretical  basis,  the  only  sites  whore  the 
instantaneous  reaction  as.sumption  may  not  apply  are  sites  with  very  low  hydiaulic  re.sidence 
times  (very  high  groundwater  velocities  iind  short  source  zone  lengths). 


K 


Distance  along  plume  centerline 


Distance  along  plume  centerline 


Figure  I.  Distribution  of  BTEX.  Electron  Acceptors,  ind  Metabolic  By-Products  vs.  Disunce  Along  Centerline 
of  Plume. 

Sampling  Dace  and  Source  of  Data:  Tyndall  3/95,  Kees/er  4/95  (Groundwater  Services,  Inc),  Patrick  3/94  (note:  I  NOS 
out/irr  removed,  sulfate  not  plotted).  Hill  7/9 J,  Elmendorf  Site  5T4I  6/94,  Elmendorf  Site  HG  10  6/94,  (Parsons 
Engineering  Suence) 
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Kinetic-limited  sites,  however,  appear  to  be  relatively  rare  as  the  instantaneous  reaction 
pattern  is  observed  even  at  sites  such  as  Siie  870  at  Hill  AFB,  with  residence  times  of  a  month 
or  less.  As  shown  in  Figure  1,  this  site  has  an  active  sulfate  reducing  and  methane  prcxluction 
zone  within  100  ft  of  the  upgradient  edge  of  plume.  With  a  1600  ft/yr  seepage  velocity  is 
considered,  this  highly  anaerobic  zone  has  an  effective  residence  time  of  23  days.  Despite  this 
very  short  residence  time,  significant  sulfate  depletion  and  methane  production  were  observed 
in  this  zone  (see  Figure  1).  If  the  anaerobic  reachons  were  significantly  constrained  by 
microbial  kinetics,  the  amount  of  sulfate  depletion  and  methane  production  would  be  much  less 
pronounced.  Therefore  this  site  supports  the  concdusion  that  the  instantaneous  reaction 
assumption  is  applicable  to  almost  all  petroleum  release  sites. 

Biodegradation  Capacity 

To  apply  am  electron-acceptor-Iimited  kinetic  model,  such  as  the  instantaneous  reaction,  the 
amount  of  biodegradation  able  to  be  supported  by  the  groundwater  that  moves  through  the 
source  zone  must  be  calculated.  The  conceptual  model  used  in  BIOSCREEN  is: 

1 .  Groundwater  upgradient  of  the  source  contains  electron  acceptors. 

2.  As  the  upgradient  groundwater  moves  through  the  source  zone,  non-aqueous  phase 
liquids  (NAPLs)  and  contammated  soil  release  dissolvable  hydrocarbons  (in  the  case  of 
petroleum  sites,  the  BTKX  compounds  benzene,  toluene,  ethylbenzene,  xylene  are 
released). 

3.  Biological  reactions  occur  until  the  available  electron  acceptors  in  groundwater  are 
consumed.  (Two  exceptions  to  this  conceptual  model  are  the  iron  reactions,  where  the 
electron  acceptor,  ferric  iron,  dissolves  from  the  aquifer  matrix;  and  the  methane 
reactions,  where  the  electron  acceptor,  CO,  is  also  produced  as  an  end-product  of  the 
reactions.  For  these  reactions,  the  metabolic  by-products,  ferrous  iron  and  methane,  can 
be  used  as  proxies  for  the  potential  amount  of  biodegradation  that  could  occur  from  the 
iron-reducing  and  inethonogenesis  reactions.) 

4.  The  total  amount  of  available  electron  acceptors  for  biological  reachons  ran  be 
estimated  by  a)  calculatmg  the  difference  between  upgradient  concentrations  and  source 
zone  concentrations  for  oxygen,  nitrate,  and  sulfate;  and  b)  mciisuring  the  produchon  of 
metabolic  by-products  (ferrous  iron  and  methane)  in  the  source  zone. 

5.  Using  .stoichiometry,  a  utilization  factor  can  be  developed  showing  the  ratio  of  the 
oxygen,  nitrate,  and  sulfate  consumed  to  the  mass  of  dissolved  hydrocarbon  degraded  in 
the  biodegradation  reactions.  Similarly,  utilization  factors  can  be  developed  to  show 
the  ratio  of  the  mass  of  metabolic  by-products  that  are  generated  to  the  mass  of 
dissolved  hydrocarbon  degraded  in  the  biodegradation  reactions.  Wiedemcier, 
Wilson,  ct  ai.  (1995)  provides  the  following  utilization  factors  based  on  the 
degradation  of  combined  BTEX  constituents: 


Electron  Acceptor/By-Product 

BTEX  Utilizrition  Factor  gm/gna 

Oxygen 

3.14 

Nitrate 

4.9 

Ferrous  Iron 

21.8 

Sulfate 

4.7 

Methane 

0.78 
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6.  For  a  given  background  concentration  ot  an  individual  electron  acceptor,  the  potential 
contaminant  mass  removal  or  "biodegradation  capacity"  depends  on  the  "utilization 
factor"  for  that  electron  acceptor.  Dividiitg  the  background  concentration  of  an  electron 
acceptor  by  its  utilization  factor  provides  an  estimate  (in  BTEX  concentration  units)  of 
the  assimilative  capacity  of  the  aquifer  by  that  mode  of  biodegradation. 

Note  that  BIOSCREEN  is  based  on  the  BTEX  utilization  provided  above.  If  other 
constituents  are  modeled,  the  utilization  factors  in  the  software  (scroll  down  from  the 
input  screen  to  find  the  utilization  factors)  should  be  changed  or  the  available  oxygen, 
nitrate,  non,  sulfate,  and  methane  data  should  be  adjusted  accordingly  to  reflect 
alternate  utilization  factors. 

When  the  available  electron  acceptor/by  product  concentrations  (No.  4)  are  divided  by 
the  appropriate  utilization  factor  (No.  5),  an  estimate  of  the  "biodegradation 
capacity"  of  the  groundwater  flowing  through  the  source  2»rie  and  plume  can  be 
developed.  ITic  biodegradation  capacity  is  tlien  used  directly  in  the  BIOSCREEN 
model  to  simulate  the  effects  of  an  instantaneous  reaction.  The  suggested  calculation 
approach  to  develop  BIOSCREEN  input  data  is: 

Biodegradation  Capacity  (mg/L)  = 

(  (Average  Upgradient  Oxygen  Cone.)  -  (Minimum  Source  Zone  Oxygen  Cone)  j  /3.14 
+  (  (Average  Upgradient  Nitrate  Cone.)  -  (Minimum  Source  Zone  Nitrate  Cone)  )  / 4.9 
+  1  (Average  Upgradient  Sulfate  Cone.)  -  (Minimum  Source  Zone  Sulfate  Cone) )  /  4.7 
+  1  Average  Observed  Ferrous  1  ■  on  Cone,  in  Source  Area}  /  21.8 

e  (  Average  Observed  Methane  Cone,  in  Source  Area  )  /  0.78 

Biodegradation  capacity  is  similar  to  "Expressed  Assimilative  Capacity"  (EAC) 
described  in  the  AFCEE  Technical  Protocol  except  that  EAC  calculations  do  not  use  the 
maximum  source  concentraHons  for  iron  and  methane.  Calculated  biodegradation 
capacities  or  EACs  at  different  U.S.  Air  Force  RNA  research  sites  have  ranged  from  7  to 
70  n-ig/l.  (.see  Table  1).  The  median  biodegradation  capacity/EAC  for  28  AFCEE  sites  is 
28.5  mg/I  . 

Note  that  one  criticism  of  tliLs  lumped  biodegradation  capacity  approach  i.s  that  it 
a-'isumcs  that  all  of  the  various  aerobic  and  anaerobic  reactions  occur  over  the  entire 
area  of  the  contaminant  plume,  and  that  the  theoretical  "zonation"  of  reactions  is  not 
simulated  in  BIOSCREEN  (e.g.  typically  dis.solved  oxygen  utilization  occurs  at  the 
downgradient  portion  and  edges  of  the  plume,  nitrate  utilization  a  little  closer  to  the 
source,  iron  reduction  in  the  middle  of  the  plume,  sulfate  reduction  near  the  .source,  and 
methane  production  in  the  heart  of  the  source  zone).  A  careful  mspection  of  acmal  field 
data  (see  Figure  1)  shows  little  or  no  evidence  of  this  theoretical  zonation  of  reactions; 
in  fact  all  of  the  reactions  appe^u•  to  occur  simultaneously  m  the  .source  zone.  The  most 
aimmon  pattern  observed  at  petroleum  release  sites  is  that  ferrous  iron  seems  to  be 
restricted  to  the  higher-concentration  or  source  zone  areas,  with  the  other  reactions 
(oxygen,  nitrate,  and  sulfate  depletion),  occurring  throughout  the  plume. 

RIOSCREHN  nssume.s  that  all  of  the  biodegradation  reactions  (aerobic  and  anaerobic) 
occur  almost  i.istantaneously  relative  to  tho  hydraulic  residence  time  in  the  .source  area 
Lind  plume,  lecause  iron  reduction  and  methane  production  appear  to  ix^iTU' only  in  the 
source  zone  (probably  due  to  the  removal  of  these  metabolic  hy-products)  it  is 
recommended  to  use  the  average  iron  and  methane  concentrations  observed  m  the  .source 
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zone  for  the  calculation  of  biodegradation  capacity  instead  of  maximum  concentrations. 

In  addition,  the  iron  and  methane  concentrations  are  used  during  a  .secondary’ 

calibration  step  (see  below).  Beta  testing  of  BIOSCREEN  Indicated  that  the  use  of  the 

maximum  concentration  of  iron  and  methane  tended  to  overpredict  biodegradation  at 

many  sites  by  assuming  these  reactions  occurred  over  the  entire  plume  area.  Use  of  an  t 

average  value  (or  some  reduced  value)  helps  match  actual  field  data. 

7.  Note  that  at  some  sites  the  instantaneous  reaction  model  wUl  appear  to  overpredict 
the  amount  of  biodegradation  that  occurs,  and  underpredict  at  others.  As  witli  the  case 
of  the  first-order  decay  model,  some  calibration  to  actual  site  conditions  is  required. 

With  the  first-order  decay,  the  decay  coefficient  is  adjusted  arbitrarily  until  the  I 

predicted  values  match  observed  field  conditions.  With  the  instantaneous  reaction 
model,  there  is  ivo  first-order  decay  coefficient  to  adjust,  so  the  following  procedure  is 
'ecommended; 

A)  The  primary  calibration  step  (if  needed)  is  to  manipulate  the  model's  ilispersivity 

values.  As  described  in  the  BIOSCREEN  Data  Entry  Section  below,  values  for  • 

dispersivity  are  related  to  aquifer  scale  (defmed  as  the  plume  length  or  distance  to 

tile  measurement  point)  and  simple  relationships  are  usually  applied  to  estimate 

dispersivities.  Gelliar  ct  al.  (1992)  cautions  that  dispersivity  values  vary  between 

2-3  orders  of  magnitude  for  a  given  scale  due  to  natural  variation  in  hydraulic 

conductivity  at  a  particular  site.  Therefore  dispersivits'  values  can  be  manipulated 

within  a  large  range  and  still  be  within  the  rimge  of  values  ob.served  at  field  test  I 

sites.  In  BIOSCREEN,  adjusting  the  transverse  dispersivity  alone  will  u.sually  be 
enough  to  calibrate  the  model. 

15)  As  a  secondary  calibration  step,  the  biodegradation  capacity  calculation  may  be 
reevaluated.  There  is  .some  judgment  involved  in  averaging  the  electron  acceptor 
concentrations  observed  in  upgradient  wells;  determining  the  minimum  oxygen,  • 

nitrate  and  sulfate  in  the  source  zone;  and  estimating  the  average  ferrous  iron  and 
methane  concentrations  in  the  source  zone.  Although  probably  not  needed  in  mo.st 
applications,  tliese  values  may  be  adjusted  as  a  final  level  of  calibration. 


« 
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BIOSCREEN  CONCEPTS 

The  BldSCKEEN  Natural  A  ienuation  software  is  ba.sed  on  the  Domenico  (1987)  three 
dimensional  analytical  solute  tran.sport  model.  The  original  model  assumes  a  fully-penetrating 
vertical  plane  source  oriented  perpendicular  to  groundwater  flow,  to  simulate  the  relea.se  of 
organics  to  movmg  groundwater.  In  addition,  the  Domenico  solution  accounts  for  the  effects  of 
advective  tran.sport,  three-dimensional  dispersion,  ad.sorption,  and  first-order  decay.  In 
BIOSCRliEN,  the  Domenico  .solution  has  been  adapted  to  provide  three  different  model  types 
representing  i)  transport  with  no  decay,  ii)  transport  with  first-order  decay,  and  iii)  tran.sport 
with  "instantaneous "  biodegradation  reaction  (see  Model  Types).  Guidelines  for  selecting  key 
Liipui  parameters  for  the  model  are  outhned  m  BIUBGKEEN  input  Baramelers.  For  help  iii 
Output,  see  13IOSCREEN  Clutput. 

BIOSCREEN  Model  Types 

The  software  allows  tile  user  to  see  results  from  three  different  types  at  groiuui  water  tran.sport 
models,  all  based  on  the  Domenico  solution; 


I 


I 
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1.  Solute  transport  with  no  decay.  This  model  is  appropriate  for  predicting  the  movement 
of  conservative  (non-degrading)  solutes  such  as  chloride.  The  only  attenuation 
mechanisms  are  dispersion  in  the  longitudinal,  transverse,  and  vertical  directions,  and 
adsorption  of  contaminants  to  the  soil  matrix. 

2.  Solute  transport  with  first-order  decay.  With  this  model,  the  solute  degradation  rate 
is  proportional  to  the  solute  concentration.  The  higher  the  concentration,  the  higher 
the  degradation  rate.  This  is  a  conventional  method  for  simulating  biodegradation  in 
dissolved  hydrocarbon  plumes.  Modeler  using  the  first-order  decay  model  typically 
use  the  hrst-order  decay  coefficient  as  a  calibration  parameter,  and  adjust  the  decay 
coefficient  until  the  model  results  match  field  data.  Witla  this  approach,  uncertainties 
in  a  number  of  parameters  (e.g.,  dispersion,  so-ption,  biodegradation)  are  lumped 
together  in  a  single  calibration  parameter. 

Literature  values  for  the  half-life  of  benzene,  a  readily  biodegradable  dissolved 
hydrocarbon,  raiige  from  10  to  730  days  while  the  half-li/e  for  TCE,  a  more  recalcitrant 
constituent,  is  10.7  months  to  4.5  years  (Howard  et  ai,  1991).  Other  applications  of  the 
first-order  decay  approach  include  radioactive  solutes  and  abiotic  hydrolysis  of 
selected  organics,  such  as  dissolved  chlorinated  solvents.  One  of  the  best  sources  of 
first-order  decay  coefficients  in  groundwater  systems  is  The  Hanabook  of 
Environmental  Degradation  Ratea  (Howard  et  ai,  1991). 

The  first-order  decay  model  does  not  account  for  site-specific  mfomration  such  as  the 
availability  of  electron  acceptors.  In  addition,  it  does  not  assume  any  biodegradation  of 
dissolved  constituents  in  the  source  zone.  In  other  words,  this  model  assumes 
biodegradation  starts  immediately  downgradient  of  the  .sourct',  and  that  it  does  not 
depress  the  concentrations  of  di-s.solved  organics  in  the  source  zone  itself. 

3.  Solute  transport  with  "instantaneous"  biodegradation  reaction.  Modeling  work 
conducted  by  GSI  indicate  first-order  expre.s.dons  may  not  be  as  accurate  for  describing 
natural  attenuation  processes  as  the  instantaneous  reaction  assumption  (Connor  et  aL, 
1994).  Biodegradation  of  organic  contaminants  in  groundwater  is  more  difficult  to 
quantify  using  a  fir.ct-order  decay  equation  because  electron  acceptor  limitations  are  not 
con.sidercd,  A  mom  accurate  prediction  of  biodegr.idation  effects  may  bo  realized  by 
incorporati.ig  the  instantaneous  reaction  equation  into  a  transport  model.  Thi  > 
appro.idi  forms  the  basis  for  the  BIOSCKEEN  instantaneous  reaction  model. 

To  incorporate  the  instantaneous  reaction  in  BIOSCREEN,  a  superposition  method  was 
used,  by  this  method,  contaminant  mass  concentrations  at  any  h  cation  and  time  within 
the  flow  field  are  corrected  by  subtracting  1  mg/I,  organic  mass  for  each  mg/I.  of 
biodegradation  capacity  provided  by  all  of  the  available  electron  acceptors,  in 
accoidance  with  the  m.stiuitaneuus  reaction  assumption.  Borden  et  al.  (1986)  concluded 
that  this  simple  superposition  technique  was  an  exact  replacement  for  more 
sophisticated  oxygen-limited  expres.sions,  as  long  as  the  oxygen  and  hydrocarbon  had 
the  same  transport  rates  (e.g.,  retardation  factor,  R  =  1).  Connor  ef  al.  (1994)  revived 
iliis  approach  for  iLse  in  spreadsheets  and  compared  the  results  to  those  from  more 
sophisticated  but  difficult  to  use  numerical  nrodels.  Tlrev  found  this  approach  to  work 
well,  even  for  retardation  factors  griMter  than  1,  so  this  superposition  approach  was 
incorporated  into  the  BK^SCREF.N  model  (see  Appendix  A. 2). 
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Which  Kinetic  Model  Should  One  Use  in  BIOSCREEN? 

BIOSCREEN  gives  the  user  three  different  models  to  chi'ose  from  to  help  see  the  eCect  of 
biodegradation.  At  almost  all  petroleum  release  sites,  biodegradation  is  present  and  can  b«; 
verified  by  demonstrating  the  consumption  of  aerobic  and  anaerobic  electron  acceptors. 
Therefore,  results  from  the  No  Biodegradation  model  are  intended  only  to  be  used  for 
comparison  purposes  and  to  demonstrate  the  effects  of  biodegradation  on  plume  migration. 

Some  key  factors  for  comparison  of  the  First-order  IDecay  model  md  the  Instanta  ieou-  Reaction 
model  are  presented  below: 


msmmm 

Able  to  Utilize  Data  from 
AFCEE  Intrinsic  Remediation 
Protocol? 

•  No  -  Does  not  account  for 

electron  acceptors/by-products 

•  Yes  -  Accounts  for  availabiliti/  of 
electron  acceptors  and  by¬ 
products 

Simple  to  Use? 

•  Yes 

•  Yr.s 

Simplification  of  Numerical 
Model? 

•  Yes  -  many  numerical  models 
tncliuie  first-order  decay 

•  Yes  -  Simplification  of 

BIOPLUME  III  model 

Familiar  to  Modelers? 

•  More  commonly  used 

•  Used  less  frequently 

Key  Calibration  Parameter 

•  First-Order  Decay  Coefficients 

*  Soiircf  Term/Dispersivitv 

Over  -  or  Underestimates 
Source  Decay  Rate? 

•  May  underpredict  rate  of 
source  depletion  (see  Newell 
etai.  1995) 

•  May  be  r-ore  accurate  for 
estimating  rate  of  source 
depletion  (see  Newell  cl  al..  1995) 

A  key  goal  of  the  AFCOE  Natural  Attenuation  Initiative  is  to  quantify  the  magnitude  of  RNA 
based  on  field  measurements  of  electron  acceptor  consumption  and  metabolic  by-product 
production.  Therefore,  the  Instantaneous  Reaction  model  is  recommended  either  alone  or  in 
addition  to  the  first-order  decay  model  (if  appropriate  calibration  is  performed)  for  mo.st  sites 
where  the  Intrinsic  Remediation  Technical  Protocol  (Wiedemeier,  Wilson,  et  al.,  1995)  has 
been  applied.  For  a  more  rigorous  analysis  of  natural  attenuation,  the  BIOPLUME  III  model  (to 
be  released  in  late  1996)  may  be  more  appropriate. 


BIOSCREEN  DATA  ENTRY 

Three  important  considerations  regarding  data  input  are: 

1  To  .see  the  example  data  set  in  the  input  screen  of  the  software,  click  on  the  "Paste 
Example  Dat<.  Set"  button  on  the  lower  right  portion  of  the  input  screen. 

2)  Because  BIOSCREHN  is  based  on  the  Excel  spreadsheet,  you  hai'c  to  click  outside  of 
the  cell  where  you  just  entered  data  or  hit  "return"  before  any  ol  the  buttons  will 
work. 

.1)  Several  cells  have  data  that  can  be  entered  directly  nr  can  be  calculated  by  the  model 
using  data  entered  in  the  grey  cells  (e.g.,  seepage  velocity  can  be  entered  directly  or 
calculated  using  hydraulic  conductivity,  gradient,  emd  effective  porosity).  If  the 
calculation  option  does  not  appear  to  work,-  check  to  make  sure  that  there  is  still  a 
formula  in  the  cell.  It  not,  you  can  restore  the  formula  hy  clicking  (n  the  "Restore 
Formulas"  huttou  in  the  bottom  right  hand  side  of  the  input  .screen.  If  there  still 
appears  to  I  a  problem,  click  somewhere  outside  of  the  last  cel!  where  you  entered 
data  and  then  click  on  tlsc  "Recalculate"  button  on  tire  input  screen. 
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JHYDRdGEOLOGIC  DATA  < 


^P*rameteri 


Typical  Values 


Source  of  Data 


Actual  interstitial  groundwater  velocity,  equalmg  Darcv  velocitv 
divided  by  effective  porosity.  Note  th.it  the  Domenico  model  and 
BIOSCREEN  are  not  formulated  to  simulate  the  effects  of 
chemical  diffusion.  Therefore,  contaminant  transport  through 
very  slow  hydrogeologic  regimes  (e.g.,  clays  and  slurrv  walls) 
should  probably  not  be  modeled  using  BIOSCREEN  unless  the 
effects  of  chemical  diffusion  are  proven  to  be  insignificant. 
Domenico  and  Schwartz  (i990)  mdicate  that  chemical  diffusion  is 
insignificant  for  Peclet  numbers  (seepage  velocitv  times  median 
pore  size  divided  by  the  bulk  diffusion  coefficient)  >  100. 


0.5  to  200  ft/yr 


Calculated  by  multiplymg  hydraulic  conductivity  by  hydraulic- 
gradient  and  dividing  by  effective  porosity.  It  is  strongly 
recommended  that  actual  site  data  be  used  .or  hydraulic 
conductivity  and  hydraulic  gradient  data  parameters;  effective 
porosity  can  be  estimated. 


1)  Enter  directly  or  2)  Fill  in  values  for  hydraulic  conductivity, 
hydraulic  gradient,  and  effective  porosity  as  described  below  and 
have  BIOSCREEN  calculate  seepage  velocity.  Note,  if  the 
calculation  option  dcxis  not  appear  to  work,  check  lo  make  sure 
that  fht  cell  still  contains  a  fomruln.  If  not,  you  can  reincarnate 
the  formula  by  clicking  cn  the  "Restore  .‘-■ormulas"  button  cn  the 
bottom  right  hand  side  of  the  input  screi  a  If  there  is  still  .i 
problem,  make  sure  to  click  somewhere  outside  of  the  last  cell 
where  you  entered  data  and  then  click  on  the  "Recalculate"  button 
on  the  input  screen. 


Parameter  ' 

H/draulic  ConducUvity ;  (K)'  .s-'--  vy'";  ’..'.liC 

Units 

cm /sec 

Description 

flon/.ontal  hvdr.iulic  conductivitv  oi  the  s-iturated  porous 
medium. 

Typical  Values 

CUiys.  <1x10''' cm/s 

Silts;  1x10'^  -  lx  11)-'*  cm/s 

Silty  sands  IxlO"'*  -  1x10-’ cm/s 

Clean  s.ind.s  IxlO’’  -  1  cm/s 

(iravel.s  >  1  cm/s 

Source  of  Data 

Pump  tests  or  slug  tests  al  the  site,  it  is  strnngly  nH  ommeniied 
that  .actual  site  data  be  used  for  most  KNA  studies 

How  to  Enter  Data 

I'.nter  directly.  If  seep.lgi-  velocity  is  enlered  directlv,  this 
p.ir.imeter  is  not  needed  in  Bit  )S( 'REPIN. 

Units 

ft/ft 

Description 

The  slope  of  the  potentiometric  surface.  In  unconfined  aquifers, 
this  is  equivalent  to  the  slope  of  the  water  table. 

Typical  Values 

0.0001  -  0.05  ft/ft 

Source  of  Data 

Calculated  by  constructing  potentiometric  surface  maps  using 
static  water  level  data  from  monitoring  wells  and  estimating  the 
slope  of  the  potentiometric  surface. 

How  to  Enter  Data 

Enter  directly.  If  seepage  velocity  is  entered  directly,  this 
parameter  is  not  needed  in  BIOSCREEN. 

Parameter  ,  “ 

Effective  Porasity 

Units 

unitless 

Description 

Dimensionless  ratio  of  the  volume  of  interconnected  voids  to  the 
bulk  volume  of  the  aquifer  matrix.  Note  that  “total  porosity"  is 
the  ratio  of  all  voids  (mc)uded  non-connected  voids)  to  the  bulk 
volume  of  the  aquifer  matrix.  Difference  between  total  and 
effective  porosity  reflect  lithologic  controls  m  pore  structure.  In 
unconsolidated  sediments  coarser  than  silt  size,  effective  porosity 
can  be  less  than  total  porosity  by  2-5'’/l>  (e.g.  0.2b  vs,  0.30)  (Smith 
and  Wheatcratt,  1993). 

Typical  Values 

Values  for  Effective  Porosity: 

Clay  0.01-0.20  Sandstone  0.005  -  0.10 

Silt  0.01  -"10  Unlract.  Limestone  0.001-  0.05 

Fine  Sand  0.10  -  0.30  Fract.  Granite  0.00005  -  0.01 

Medium  Sand  0.15-l)..30 

Coarse  Sand  0.20  -  0.35 

Gravel  U.  10 -0.35 

(Frimi  VV,.-ilcmt’itT,  Wilson.  Il'rom  Domenico  anJ  Schwariz.  IIIHOl 

L't  dl.,  >'9Vt),  anyiriu//v  /rom 

Domenico  anJ  S-  hwartz,  1990 
and  Walton.  I'lHtt) 

Source  of  Data 

Typically  estimated.  (Sie  commonly  used  value  tor  silts  and  sands 

IS  an  effective  porosity  of  0.25.  Tlie  ASTM  RBCA  Standard 
(ASTM,  1995)  includes  a  default  value  of  ().3H  (to  fx-  lekd 
pf!marily  for  unconsolidated  deposits). 

I'l’ittF  dircCilV-  iMrti  if  VOiuCiiV  iS  CntOifii  uiiL'Lily, 

this  parameter  is  still  needed  to  i:ak;ulate  the  retardation  fr.rtor 

and  plumt?  mass. 
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2.  PISPCRISIVITY 


Units 

tt 

Description 

Dispersion  refers  to  the  process  whereby  a  plume  will  spread  out  in  a 
lon^titudinal  direction  (alonj;  the  direction  ol  groundwater  flow), 
transversely  (perpendicular  to  groiindwater  tlow),  and  vertically 
downivards  due  to  mechanical  mixing  in  the  aquifer  and  chemical 
dituisinn  Si’loction  ot  dispersivitv  values  is  a  difficult  process, 
gnen  the  mipracticahihtv  ot  measiinng  dispeisum  in  the  field 
Howcmt.  simple  estimation  techniques  based  m  the  length  ot  the 
plume  or  distance  to  the  measurement  point  ("scale")  are  available 
irom  a  compilation  ot  field  test  data  Note  that  researchers  indicate 
that  dispersivitv  values  can  range  over  2-T  orders  ol  magnitude  tor  a 
given  value  ig  plume  length  or  distance  to  measurement  point 

K  lelhar  rt  ii/  .  W02',  In  HKfSCRIihN,  dispersu  itv  is  uvxf  as  the 
primary  calibration  parameter  (see  pg  U').  h'or  mon'  intorniLition  in 
dispersivitv.  see  Appendix  A.. 3,  pg  40). 

Typical  Values 

Typical  dispersivitv  relationships  as  a  tunction  ot  1  p  (plume  length 
or  distance  to  measurement  point  in  It)  are  provided  below 
HUWC'KT'h'N  IS  programmed  with  some  ci>mmonlv  ustxi  relationships 
representative  ot  tvpical  and  low-end  dispersivities 

•  Longitudinal  LTispersivity 

Alphas  P2X  O.Xt  t.V.i  .111./  IcDt, 111.  1  as,, 

(1  ,  III  til 

•  Transverse  Dispersivity 

•Alph.l  V  ■  0.10  .llpha  s  i/liiv-./  an  /iig/i  niiiiliihJy 

/ii'iiils  tiiiiii  1  .liliur  I’t  ill  .  I'l')  ’ ' 

•  Vertical  Dispersivitv 

.Alph.l  /  cerv  low  (i.e.  1  \  10  "Oj  i/t.oo/ an  .  lUisci.-.Kiff  ,--,linuitr‘ 

(  •  fifniTinnK  um’M  n*l.iUf'ir%hn'-*  ii\i  luilr 

Alpli.f  X  1)  1  1  [>  ‘ /'ji  i  tM  mA. 

Aljili.i  V  (}  l\-\. 

Alph.l/  0.(Jt  .ilph.i  \  (A'.ilXy 

Alph.l  /  0  Ii2  *  .il|’li.i  \  tt  1 1)  1  .ilph.i  V  '  M'. \ 

Source  of  Data 

Ivpic.lllv  aslim.lted  ilsing  the  re l.il uii \  .1  up-,  piovuled  ,ilio\  e  (see 
,\('pandis  .A  !,  [)g  (0) 

How  to  Enter 

Data 

1  1  T.nter  direi  tiv  oi  J|  Till  iii  v.ihie  ot  the  estim.lled  plume  leng.lh  .iiui 
h.ive  hit  )'■>( 'KT  I'N  1. ill  111. lie  till'  dispai  su  llies 

>«• 


n 


•  • 


» 


» 


r; 


» 


BlOSCREEN  U'ipr’s  Manual 


!une  1996 


Units 

It 

Description 

Estimated  length  (in  feet)  of  the  existing  or  hypothetical 
groundwater  plume  being  modeled.  This  is  a  key  parameter  as  it  is 
generally  used  to  estimate  the  dispersivity  terms  (dispersivity  is 
difficult  to  measure  and  field  data  are  rarely  collected). 

Typical  Values 

For  B'FEX  plumes,  50  -  500  ft.  For  chlorinated  solvents,  50  to  1000  ft. 

Source  of  Data 

To  simulate  an  actual  plume  length  or  calibrate  to  actual  plume  data, 
enter  the  actual  length  of  the  plume.  If  trying  to  predict  the  maximum 
extent  of  plume  migration,  use  one  of  the  two  methods  below, 

1)  Use  seepage  velocity,  retardation  factor,  and  simulation  time  to 
estimate  plume  length.  While  this  may  underestimate  the  plume 
length  for  a  non-degrading  solute,  it  may  overestimate  the  plume 
length  for  either  the  first-order  decay  model  or  instantaneous  reaction 
model  if  biodegradation  is  significant. 

2)  Estimate  a  plume  length,  nm  the  model,  determine  how  long  the 
plume  is  predicted  to  become  (this  will  vary  depending  cn  the  type  of 
kmetic  expression  that  is  u.sed),  reenter  this  value,  iind  then  mrun  the 
model.  Note  that  con'-iderable  time  and  effort  can  be  expended  trying 
to  ad)Ufit  the  estimated  plume  length  term  to  match  exactly  the 
predicted  modeling  length.  In  practice,  most  modelers  make  the 
assumption  that  dispersivity  values  are  not  very  precise,  and 
therciore  select  ball-park  values  based  in  estimated  plume  lengths 
that  are  probably  i  25"o  of  the  actual  plume  length  used  in  the 
simulations.  Note  that  BlOSCREEN  is  very  sensitive  to  the 
dispersion  estimates,  particularly  tor  the  instantaneous  reaction 
model. 

How  to  Enter 

Finter  directlv.  It  dispersivity  data  are  entered  directly,  this 

Data 

p.irameter  is  not  needed  in  HIOSC’REEN. 

3.  ADSORPTION  PATA 


Units 


Description 


Typical  Values 


Source  of  Data 


The  rate  at  which  dissolved  contaminants  movmg  througli  an 
aLjuitor  can  be  teduccd  by  sorption  of  contaminants  to  the  solid 
aquifer  matris.  Fhc  degree  of  retardation  depends  on  both  aquifer 
and  constituent  properties.  The  retardation  factor  is  the  ratio  of 
the  groundwater  seepage  velocity  to  the  rate  that  organic 
chemicals  migrate  in  the  groundwater.  A  retardation  value  of  2 
indicates  that  it  the  groundwater  seepage  velocity  is  UX)  ft/yr, 
then  the  organic  chemicals  migrate  at  approximately  50  ft/yr. 

BIOSCKHHN  simulations  using  the  instantaneou.s  reaction 
assumption  a‘  sites  with  retardation  factors  greater  than  n  .should 
be  performed  with  caution  and  verified  asing  a  moiv 
sophisticated  model  such  as  BIOlT.UMli  III  (see  Appendix  A. 2). 


1  to  2  (tor  BTI'X  in  typical  shallow  aquitci  ,) 


Usiiallv  estimated  from  soil  and  chemical  data  vrsing  v'ariablcs 
described  below  (ph  -  bulk  density,  n  ■  porosity,  K.oc  ^  organic 
carbon-watcr  partition  coefficient,  Kd  =  distribution  coefficient, 
and  toe  =  traction  org.iiuc  carbon  on  uncont.iniinated  soil)  with  the 
tollowing  expression: 


R  ■  \  \ 


K.l  ■  ph 


How  to  Enter  Data 


1m  sjiinc  cases,  the  lel.irdation  tacloi  can  tie  estimated  by 
(lanparing  the  length  ol  a  plume  afiected  by  adsorption  (such  as 
ihc  bcn/cnc  plume)  with  the  length  ot  plume  that  is  not  affected 
bv  adsorption  (such  as  chloride).  Most  plumes  do  not  have  both 
typos  ot  contaminants,  vi  it  is  mon'  aunmon  to  iisi-  the  estimation 
teclinu|ue  (see  data  entry  boxes  below). 


1)  I'litei  directly  or  2)  I'll  in  the  estimated  value.s  lor  bulk 
densilv,  partition  eoi-ttu  l•'nt,  and  tr.iclion  organic  cartion  as 
ilcscnht'd  below  and  have  liK  )S( 'Klil-N  c.ilculate  retardation. 


Parameter  ■ 

Soil  Bulk  Density  (p  i,) 

Units 

kg/I.org/niT 

Description 

bulk  densitv,  in  kg/i  .  ot  the  .iquih  r  niatrix  (related  to  porosity 
and  pure  snluls  density). 

Typical  Values 

Although  this  value  can  Ih-  uieasiiri  1  in  the  lab,  in  most  cases 
estimated  values  ,ire  usetl  A  value  ot  1.7  kg/I,  is  used  treqiieiltlv. 

Source  of  Data 

l  itliei  from  an  .in.ilvsis  nl  sod  sample-,  at  a  geoteehnu  al  lab  or 
more  eommonlv,  application  ot  estimated  values  such  a-i  i.7  kg/I.. 

How  to  Enter  Data 

I  nter  dirci  tlv  11  the  let.ird.itioii  lacUii  is  enteied  directlv,  this 
p.irameter  is  not  nei-dixi  in  BIOS(  KliliN. 

Units 

(mg/ kg)  /  (mg/L)  or  (L/kg)  or  (mL/ g) 

Description 

Chemical-spec  fic  partition  coefficient  between  soil  organic  cartion 
and  tire  aqueous  phase.  Larger  values  indicate  greater  affinity  of 
contaminants  for  the  organic  carbon  fraction  of  soil.  This  value  is 
chemical  specific  and  can  be  found  in  chemical  reference  books. 
Note  that  many  usere  of  BIOSCREEN  will  simulate  BTEX  as  a 
single  coiistituent.  In  this  case,  either  an  .average  value  for  the 
BTEX  compounds  can  be  used,  or  it  can  be  assumea  that  all  of  the 
BTEX  compounds  have  the  same  mobility  as  benzene  (the 
constituent  with  the  highest  potential  risk  to  human  health). 

Typical  Values 

Benzene  38  L/kg  Lthylbenzene  95  L/kg 

Toluene  135  L/kg  Xvlene  240  L/kg 

(ASTM,  1995) 

(Note  that  there  is  a  wide  range  of  reported  values;  for  example, 
Mercer  and  Cohen  (1990)  report  a  Koe  for  benzene  of  83  L/kg 

Source  of  Data 

Chemical  reference  literature  or  relationships  between  Koe  and 
solubility  or  Koe  and  the  octanol-wator  partition  eoeffieient 
(Kow). 

How  to  Enter  Data 

Enter  directly.  If  the  retardation  factor  is  entered  directly,  this 
parameter  is  not  needed  in  BIOSCREEN. 

Parameter 


Units 


Description 


Typical  Values 


Source  of  Data 


How  to  Enter  Data 


Fraction  ^Organic  Carbon  (foe)  ■ 


unitiobs 


l-raction  of  tho  aquifer  soil  matrix  comprised  of  natural  organic- 
carbon  in  uncontaminated  areas.  More  natural  organic  carbon  means 
higher  ad.sorption  of  organic  constituents  on  the  aquifer  matrix. 


0.0002  -  0.02 


I  he  traction  organic  carbon  value  should  be  measured  if  possible  bv 
collectmg  a  sample  of  acjuiler  inatenal  Irom  an  uneontaminated 
/.one  and  pertonning  a  laboratory  analysis  (e.g.  ASTM  Method 
2474-87  or  equivalent).  If  unknown,  a  detault  value  of  U.tXfl  e.  otteri 
used  (e.g.,  ASTM  1445). 

b.iiter  dir"ctlv  It  thir  retanlation  factor  is  entered  directly,  this 
parameter  is  not  iv'edctl  in  lllOSCRlihN. 


Typical  Values 


Rato  coefficient  describing  first-order  decay  process  for  dissolved 
constituents.  The  first-order  decay  coefficient  equals  0.693  divided 
by  the  half-life  of  the  contaminant  in  groundwater.  In 
BIOSCREEN,  the  first-order  decay  process  as.sumes  that  the  rate 
of  biodegradation  depends  only  cn  the  concentration  of  the 
contaminant  and  the  rate  coefficient.  For  example,  consider  3  mg/I. 
benzene  dissolved  m  water  in  a  beaker.  It  the  half-life  of  the 
benzene  in  the  beaker  is  72b  davs,  then  the  concentration  of  benzene 
728  davs  from  now  will  be  1.5  .iig/L  (ignonng  volatilization  and 
other  loss.  ). 

Considerable  care  must  be  extrci.sed  m  the  selection  of  a  first-order 
decav  coefficient  for  each,  constituent  in  order  to  avoid 
significantly  over-predicting  or  under-predicting  actual  decay 
rates.  Note  that  the  amount  of  degradation  that  rxeurs  is  related 
to  the  time  the  contaminants  spend  m  the  aquifer,  and  that  this 
parameter  is  not  related  to  the  time  it  takes  for  the  soume 
concentrations  to  decay  by  half. 


0  1  to  3b  yr  ‘  (see  half-life  values) 


Optional  methods  for  selection  ot  appropriate  decay  coefficients 
are  as  follows: 

Literature  Values;  Various  published  references  are  available 
listing  decay  half-life  values  for  hvdrolysis  and  biodegradation 
(e.g.,  see  Howard  ft  «/.,  1991).  Note  that  many  references  report 
the  half-lives;  tfiese  values  can  tx'  converted  to  the  fir.st-order 
decay  coefficients  using  k  -  t).693  /  t,  .  (see  dissolved  plume  half- 
life). 

Calibrate  to  Existing  Plume  Data;  It  the  plume  is  in  a  steadv-state 
or  diminishing  condition,  BKiSC'KFFN  can  Ise  u.ms.1  to  determine 
first order  decav  coefficients  that  best  m.itch  the  observed  site 
concentration-,  (hie  may  adopt  a  tri.il  .iiid-error  pr(X’t‘dure  to 
derive  a  best-fit  decav  coefficient  for  each  contaminant  For  still- 
exp.inding  plumes,  this  steailv-state  calibr.ilion  method  m.iv  over¬ 
estimate  actual  decav  rate  coetticients  and  contribute  to  an  under¬ 
estimation  of  preiiicted  concentration  levels 


11  l.nter  directly  or  ?)  Fill  in  tfie  estini.iteil  half-Ilf'.'  '.allies  as 
desciilifd  below  and  have  BK'SF  KI  FN  calculate  the  first-order 
decav  'oeflicieilts 


How  to  Enter  Data 


*»r»nieter 


)rwpt!j:i<tii?lumg}5olut8H« 


Description 


Time,  in  years,  for  dissolved  plume  concentrations  to  decay  by  '  ie 
half  as  contaminants  migrate  through  the  aquifer.  Note  that  wie 
amount  of  degradation  that  occurs  is  related  to  the  time  the 
contaminants  spend  in  the  aquifer,  and  that  the  degradation  IS 
NOT  related  to  the  time  it  takes  for  the  source  concentrations  to 
decay  by  half. 


Modelers  using  the  first-order  decay  model  typically  use  the  first- 
order  decay  coefficient  as  a  calibration  parameter,  and  adjust  the 
decay  coefficient  until  the  model  results  match  field  data.  With 
this  approach,  uncertainty  in  a  number  of  parameters  (e.g., 
dispersion,  sorption,  biodegradation)  are  lumped  together  in  a 
single  calibration  parameter 

Considerable  care  must  be  exercised  in  the  selection  of  a  first-order 
decay  coefficient  for  each  contammant  in  order  to  avoid 
significaiitlv  over-predicting  or  under-predicting  actual  decay 
rates. 

Typical  Values  Benzene  0.02  to  2.0  yrs 

Toluem  0.02  twO.r/yr 

Ethylbenzene  U.uln  to  0.62  yr 
Xylene  0.C38  to  1  yr 

(from  ASTM.  1995) 

Source  of  Data  Optional  methods  for  selection  of  appropriate  decay  coefficients 

ere  as  follows: 

Literature  Values:  Various  published  references  are  available 
listing  decay  half-life  values  for  hydrolysis  and  biodegradation 
(e  g.,  see  Howard  cl  i  '  1991). 

Calibrate  to  Existing  Plume  Data:  If  the  plume  is  in  a  steady-state 
or  diminishing  condition,  BIOSCREEN  con  be  used  to  determine 
first-order  decay  coefficients  that  best  match  the  observed  site 
concentrations.  A  tnal-and-c  ror  prixedure  may  be  adopted  t(' 
derive  a  best  fit  decay  coellicient  for  each  contaminant.  I'or 
expanding  plumes,  this  steadv'-state  calibration  method  may  over¬ 
estimate  actual  decay-rate  coefficients  and  contribute  to  an  under¬ 
estimation  of  predicted  concentration  levels. 


How  to  Enter  Data 


Enter  directly  If  the  first-order  decay  coefficient  is  entered 
directly,  this  parameter  is  not  needed  in  BIOSCREHN 


Units 

irg/'L 

Description 

This  parameter,  used  in  the  instantaneous  reaction  model,  is  otk? 
component  of  the  total  biodegradation  capacity  of  the  groundwater 
as  It  flows  through  the  source  zone  and  contaminant  plume.  The 
model  assumes  that  3.14  mg  of  oxygen  are  required  to  ninMime  1  mg 
of  BTEX  (Wiedeineier,  Wilson,  et  al.,  1995).  Note  that  this 
parameter  is  used  for  the  instantaneous  reaction  model,  which  is 
appropriate  only  for  readily  biodegradable  compiounds  such  as 
BTEX  that  degrade  according  to  the  assumed  BIOSCREEN 
utilization  factors,  and  is  not  appropriate  for  more  recalcitrant 
compounds  such  as  the  chlorinated  solvents 

Typical  Values 

Data  trom  28  AFCEE  sites  (see  Table  1): 

Median  =  5.8  mg/I.  Maximum  =  12.7  mg/1.  Minimum  =  0.4  mg /I, 

Source  of  Data 

For  planning  studies,  typical  values  taken  from  Table  1  can  be  used. 
Fcr  actual  RNA  studies,  the  Air  Force  Intrinsu  Remediation 
Technical  Protocol  (Wiedemeier,  Wilson,  et  al  ,  1995)  should  be 
applied.  Enter  the  average  background  concentration  of  oxygen 
minus  the  lowe.st  observed  concentration  of  oxygen  in  the  source 
area.  BlOSCiREEN  automatically  applies  the  utilization  factor 
used  to  compute  a  biodegradation  capacity. 

How  to  Enter  Data 

Enter  directly. 

- 1 

Units 

mg/L 

Description 

This  parimeter,  used  m  the  instantaneous  reaction  model,  is  one 
component  of  the  total  biodegradation  capaciEy  of  the  grouiidwater 
as  It  flows  through  the  source  zone  and  contamuiant  plume.  The 
motfel  assumes  that  4.9  mg  of  nitrate  are  required  to  consume  1  mg  of 
BTF.X  (Wiedemeier,  Wilson,  et  at.,  1995).  Note  that  this 
parameter  is  used  for  th('  instantaneous  reaction  model,  which  is 
appropri.ite  only  for  readily  biodegr.id.ible  conifxiunds  .sui'h  as 
BTEX  that  degradi'  accordmg  to  the  .is.suMud  BlkiSCKEEN 
utilization  factors,  and  's  not  appropriate  for  more  recalcitrant 
compounds  such  as  the  chlorirated  soivenls. 

Typical  Values 

Data  from  28  AFCEf-.  sites  (see  1  able  1): 

Median  -  6.3  mg/L  Maxuiium  =  69.7  mg/I.  Minimum  -  0  mg/I. 

Source  of  1  ta 

For  pldiuiing  studies,  typical  values  taken  from  Table  1  can  be  uwd. 
For  actual  RNA  studies,  the  Air  Force  Intrinsic  Remediation 
Technical  ProtiKol  (Wiedemeier.  Wilson,  t'f  ul.,  1995)  should  N 
applied  Enter  the  average  backgrocmd  concentration  of  nitrate 
minus  the  lowest  observed  concentration  of  nitrate  in  lire  .siunv 
.irea  BIOSCKEEN  automatically  applies  the  utilization  faclor  to 
'impute  a  biodegradation  capacity. 

din'ctiy. 


How  to  Enter  Data 
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Units 

n-ig/L 

Description 

This  parameter,  used  in  the  instantaneous  reaction  model,  is  one 
component  of  the  total  biodegradation  capacity  of  the  groundwater 
as  it  flows  through  the  source  zone  and  contaminant  plume.  Ferrous 
iron  is  a  metabolic  by-product  of  the  anaerobic  reachon  where  solid 
ferric  iron  is  used  as  an  electron  acceptor.  The  model  assumes  that 
21.8  mg  of  ferrous  iron  represents  the  consumption  of  1  mg  of  BTEX 
{Wiedemeicr,  Wilson,  et  al.,  1995).  Note  that  this  parameter  is 
used  for  the  instantaneous  reaction  model,  which  is  appropriate 
only  for  readily  biodegradable  compounds  such  as  BTEX  that 
degrade  according  to  the  assumed  BIOSCREEN  utilization  factors, 
and  is  not  appropriate  for  more  recalcitrant  compounds  such  as  the 
chlorinated  solvents. 

Some  researchers  suggest  the  observed  ferrous  iron  concentration  is 
much  less  (107o  or  less)  tlian  tlie  actual  amount  of  ferrous  iron  that 
has  been  generated  due  to  the  sorption  of  ferrous  iron  onto  the 
aquifer  matrix  (Lovely,  1995).  If  this  is  the  case,  then  the  value 
used  for  this  parameter  should  be  mucli  higher  than  the  observed 
maximum  concentration  of  ferrous  iron  in  the  aquifer. 

Typical  Values 

Data  from  28  AFCEE  sites  (see  Table  1): 

Median  -  I6.bmg/I.  Maximum  =  599.5 mg/ L  Minimum  =  0 mg/L 

Source  of  Data 

For  planning  studies,  typical  values  taken  from  Table  1  can  be  used. 
For  actual  RNA  studies,  the  Air  Force  Intrinsic  Remediation 
Technical  Protocol  (Wiedemeicr,  Wilson,  et  al.,  1995)  should  be 
applied.  Enter  the  average  observed  concentration,  in  mg/L,  of 
ferrous  (dissolved)  iron  found  in  the  .si'urce  area  (approximately 
til,’  area  where  ferrous  iron  has  been  observed  in  monitoring  wells). 
BIOSCREEN  automatically  applies  the  utilization  factor  to 
compute  a  biodegradation  capacity. 

How  to  Enter  Data 

Enter  directly. 

.;i.  i 

livaiE-oateitaai . LI  i'  dii  / 

Units 

mg/L 

Description 

This  parameter,  used  in  the  instantaneou;.  reaction  model,  is  one 
component  of  the  total  biodegradation  capacity  of  the  groundwater 
as  it  flows  through  the  .'siurce  zone  and  contaminant  plume.  The 
model  assumes  that  1.7  mg  of  sulfate  are  rei.[uired  to  consume  1  mg  of 
BTEX  (Wiedemeier,  Wilson,  et  ul.,  1995).  Note  that  this 
parameter  is  used  -or  the  mstantaneous  reaction  model,  which  is 
appropriate  only  for  readily  biodegradable  compiunds  sudi  as 
BTEX  that  degrade  according  to  the  assumed  BIOSCRF.HN 
utilization  factors,  and  is  nc‘  appiopriate  for  moie  recalcitrant 
compounds  such  as  the  chlorinated  solvents. 

Typical  Values 

Data  from  28  AFCEE  sites  (see  Table  1): 

Median  24.6  mg/L  ■  Maximum  -  109,2  mg/L  Minimum  -  0  mg/L 

Source  of  Data 

Fur  planning  studies,  typical  values  taken  from  Table  1  can  be  used. 
For  actual  RNA  studies,  the  Air  Force  Intrinsic  Remediation 
Technical  Protocol  (Wiedemeier,  Wilson,  ft  al.,  1995)  should  be 
applied.  Enter  the  average  background  concentration  of  sulfate 
minus  the  owest  observed  concentration  of  sulfate  in  the  source 
area.  BIOSCREEN  then  computes  a  biodegradation  capacity. 

How  to  Enter  Data 

Enter  directly. 

Units 


Typical  Values 


m.y  » 


This  parameter,  ustKi  in  the  instantaneous  reaction  model,  is  one 
component  of  the  total  biodegradation  capacity  of  the  groundwater 
as  It  flows  through  the  sour'e  zone  and  contaminant  plume. 
Mi'thane  is  a  metabolic  by-product  of  rnethanogemc  activity.  The 
model  assumes  that  0.78  mg  of  methane  represents  the  consumption 
ot  I  mg  of  BTEX  (Wiedtmeier,  Wilson,  t’f  ill.,  1995).  Note  that  this 
parameter  is  used  for  the  instantaneous  reaction  model,  which  is 
appropriate  only  for  readilv  biodegradable  compounds  such  as 
FiTIiX  that  degrade  according  to  the  assumed  HIOSCREEN 
utilization  f.ictors,  and  is  not  appropriate  for  more  recalritrant 
compounds  such  as  the  chlorinated  .solvents. 


Data  from  28  AFCEH  sites  (see  Table  1): 

Median  =  7.2mg/L  Maximum  =  48  4  mg/ L  Minimum  =  0.0  mg/ L 


1  or  plannirij^  studiGD/  typical 


taken  from  Table  1  can  be  used. 


For  actual  RNA  studies,  the  Air  Force  Intrinsic  Remediation 
Technical  Protocol  (Wiedemeier,  Wilson,  ft  al.,  1995)  .should  be 
applied.  Enter  the  average  c  bserved  concentration  ot  methane 
found  in  the  source  area  (appro  cimately  the  area  where  methane 
is  observed  in  monitoring  Wi  lls).  HI(9SC'KHHN  automatically 
computes  a  biodegradation  capacity. 


Enter  directly 


How  to  Enter  Data 


L59£ 


5.  GENERAL  DATA 


,v:7- 


Description 


Typical  Values 


Source  of  Data 


How  to  Enter  Data 


Physical  dimensions  (in  feet)  of  the  rectangular  area  to  be 
moiielcd.  To  determine  contaminant  concentrations  at  a  particular 
point  along  the  centerline  of  the  plume  (a  axnmori  approach  for 
most  risk  assessments),  enter  this  distance  in  the  "MoHoled  Area 
Length'  box  and  see  the  results  by  clicking  in  the  "Run  Centerline" 
button. 

If  one  is  interested  in  more  accurate  mass  calculations,  make  sure 
most  of  the  plume  is  within  the  zone  delineated  by  the  Modeled 
Area  Length  and  Width,  l-'ind  the  mass  balance  result.s  using  the 
"Run  Arrav"  button. 


10  to  1000  ft 


Values  should  be  slightly  larger  tlian  the  final  plume  dimensions 
or  .should  extend  to  the  downgradient  point  of  concern  (e,g.,  point  of 
exposure).  If  only  the  centerline  output  is  used,  the  plume  width 
parameter  ha.s  no  effect  on  the  results. 


Enter  dircctiv. 


Parameter 


Simulation  .Time  (t)  ■  s. 


Description 


Typical  Values 


Source  of  Data 


Time  (m  years)  for  which  concentrations  are  to  be  calculated.  E'or 
steady-state  simulations,  enter  a  large  value  (i.e.,  1CHX1  years 
would  be  sufficient  for  most  sites). 


1  to  1000  years 


To  match  an  existing  plume,  estimate  the  time  between  the 
original  release  and  the  date  the  field  data  were  colloeted.  To 
predict  the  niaxiiTuim  extent  of  plume  migration,  increase  the 
smuilalion  time  until  the  plume  no  longer  increases  in  lengt.h. 


How  to  Enter  Data 


Enter  directly. 


Units 


Description 


The  r>crr.cnii;L>  (l^^S?)  model  assumes  a  vertical  plane  source  of 
constant  concentration.  For  many  fuel  spill  sites  the  thickness  of 
this  soun-c  zone  is  only  5-20  ft,  as  petroleum  fuels  ari'  LNAPl  s 
(light  non-atiueous  phase  liquids)  that  float  on  the  water  table. 
Therefore,  the  residual  source  zones  that  are  slov/lv  dissolving, 
creating  the  dissolved  BTFiX  plume,  are  typically  restricted  to  the 
upper  part  of  the  aquifer. 


Surface 


Top  of  Water- 
Bearing  Unit 

Source  Thickness 


Bottom  of  Water- 
Bearing  Unit 


Typical  Values 

5-50  ft 

Source  of  Data 

This  value  is  usually  detennined  by  evaluating  groundwati-r  data 
from  wells  near  tlie  source  zone  screened  at  different  depths.  It  this 
type  of  information  is  not  available,  then  one  could  estimate  the 
amount  of  water  table  fluctuation  that  has  occurred  since  the  time  of 
the  relea.se  and  use  this  value  as  the  source  zone  thickness  (eijiialing 
to  the  smear  zone).  Otherwi-.e,  a  simple  assumption  of  10  feet  would 
probably  be  appropriate  tor  many  petroleum  release  sites.  Note 
that  if  DNAPl.s  are  pre.senl  at  the  site  |e.g.,  a  chlorinated  solvent 
site),  a  larger  source  zone  thickness  would  probably  be  required 

How  to  Enter  Data 


Description  The  Domenico  (1987)  mode!  assumes  a  vertical  plane  source  of 

constant  concentration.  BIOSCREEN  expands  the  simple  one  source- 
zone  approach  by  allowing,  up  to  five  source  zones  with  different 
concentrations  to  account  for  spatial  variations  in  the  source  area. 

Typical  Values  10  -  200  ft 

Source  of  Dat.  To  define  a  varying  source  concentration  across  the  site: 

1)  Draw  a  line  perpendicular  to  the  grouitdwater  flow  direction  in 
the  source  zone.  The  souire  zone  is  typically  defined  as  being  the 
area  with  contaminated  soils  having  high  concentrations  of  sorbed 
organics,  free-phase  NAPLs,  or  residual  NAPLs.  If  the  source  zone 
covers  a  large  area,  it  is  best  to  choose  the  most  downgradicnf  or 
widest  point  in  the  source  area  to  draw  the  perpendicul  ir-to-flow 


2)  Divide  the  line  into  1,  3,  or  5  zones.  A  total  of  5  zones  is  shown  on 
the  input  screen. 

3)  Determine  the  width  and  corresponding  average  concentration  ot 
Zones  1,  2,  and  3.  Typically  Zone  3  will  contain  the  highest 
concentration.  Note  that  the  model  a.ssumes  the  source  zone  us 
symmetrical  and  will  automatically  defme  source  zones  4  and  5  to  be 
identical  to  Zones  2  and  1.  Therefore,  it  i.s  not  necessary  to  specify  a  11 
.“i  zones.  For  simpler  problems,  you  can  either  use  three  zones  to  define 
varying  source  concentrations  across  the  site  (enter  information  in 
Zones  2  mid  3,  and  the  model  will  define  Zone  4)  or  just  use  a  single 
zone  (enter  data  to.-  Zone  3  only). 

4)  Enter  the  vvidth  and  source  concentration  into  the  appropriate 
zones  (>n  the  spreadsheet  For  example,  if  a  total  source  width  of  100 
ft.  is  divided  into  five  zones,  enter  20  ft  for  each  z,one  width.  Enter 
the  average  concentration  observed  acro.is  each  zone. 


Surface 


Top  of  Water- 
Bearing  Unit 


^  Need  Width  and 

I  Concentration 

I  /  of  Source  Zones 


\  \ 


Bottom  of  Water- 
Bearing  Unit 

How  to  Enter  Data  l-nter  direetlv. 
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mgA- 

Description 

BIOSCREEN  requires  source  zone  concentrations  that  correspond  to 
the  source  zone  width  data  (see  previous  page).  Suggested  rules  of 
thumb  regarding  how  to  handle  multiple  constituents  are: 

1) If  the  maximum  plume  length  is  desired,  model  liunpied 
constituents  (such  as  BTE.\).  If  a  risk  assessment  is  being  performed, 
data  on  individual  constituents  arc  needed. 

2)  If  lumped  constituents  are  being  mo.  N'led  (BTEX  all  together), 
use  either  average  values  for  the  chemical-specific  data  (Ki'c  and 
lambda)  or  the  worst-case  values  (e.g.,  use  the  lowest  of  the  Koc  and 
lambda  from  the  group  of  constituents  being  modeled)  to 
overestimate  coiiLenlrations.  Most  modeling  will  be  performed 
assuming  that  the  ratio  of  BTEX  at  the  edge  of  the  plume  is  the- 
same  as  at  the  source.  For  more  detailed  modeling  studies,  Wilson 
(1996)  has  propo.sed  the  following  rules  o  help  account  for  difterent 
rates  of  reaction  among  the  BTEX  compounds: 

•  It  the  site  is  dominated  by  aerobic  degradation  (most  i>f  the 

biodegradation  capacity  is  from  oxygen,  a  relatively  rare 

cxrairroncc)  assume  that  the  benzene  will  degrade  first  -and  that 
the  di.ssolved  material  at  the  edge  ot  the  plume  i.s  primarily  TEX. 

•  It  the  site  is  dominated  by  nitrate  utilization  (most  ot  the 

biodegradation  capacity  is  from  nitrate,  a  relatively  rare 

occuiTcnce)  assume  that  benzene  will  degrade  last  and  that  the 
dissolved  niaierial  at  the  edge  of  the  plume  is  primarily  benzene. 

•  It  the  site  is  domuiated  by  sulfate  reduction  (most  ot  the 

biodegradation  c.apacity  i.s  due  sulfate  utilization,  a  more 

common vxxunence)  assume  that  the  fx-nzene  will  degrade  at  the 
same  rate  as  the  TEX  constituents  and  that  the  dissolved  material 
at  the  edge  ot  the  plume  i.s  a  mixture  ol  BTE.X. 

•  It  th<'  site  IS  dominated  bv  metliane  production  (most  ot  the 
biodegradation  capacity  i.s  due  to  methanogene.sis,  a  mun'  c  inmoii 
ixcurrence)  assume  that  benzene  will  degrade  last  and  tfi  t  the 
dissolved  material  at  the  edge  of  the  plume  is  primarily  benzene. 

,i)  If  individual  constituents  .ire  being  modeled  with  the 
instantaneous  re.ietion  assumption,  note  that  the  total 
biodegradation  capacity  mast  lx-  reduced  to  account  for  electron 
acceptor  utilization  by  other  constituents  present  in  the  plume,  l  or 
example,  in  order  to  model  benzene  as  an  individual  constituent 
using  the  instantaneous  reaction  model  in  a  BTEX  plume  cont.iining 
equal  source  concentrations  of  benzene,  tolueme,  ethylbenzene  and 
xylene,  the  amount  ol  oxygen,  nitrate,  sulfate,  iron,  and  metliane 
should  be  reduced  by  75'/f.  to  account  for  utiliz.ition  by  tolui-ne, 
ethylbenzene,  and  xylene. 

Typical  Values 

().ni0tol20mg/I. 

Source  of  Data 

Source  area  monitoring  well  data  (see  figure  on  previous  p.ige). 

How  to  Enter  Data 

Enter  directly. 

Description  The  Domenico  (1987)  model  assumes  the  source  is  infinite,  i.e.  the  source 
concentrations  are  constant.  In  BIOSCREEN,  however,  an  approximation  for 
a  declining  source  concentration  has  been  added.  Note  that  this  is  an 
experimental  relationship,  and  it  should  be  applied  with  caution.  The 
declining  source  tcnn  Ls  ba.sed  on  the  following  assumptions: 

‘  There  is  a  finite  mass  ot  organics  in  the  source  /.one  present  as  a  free-phase 
or  residual  NAPL.  lire  NAl’L  in  the  source  zone  dissolves  slowly  as  fresh 
groundwater  passes  through. 

•  The  change  in  source  /.one  concentration  cait  be  approximated  as  a  first- 
order  decay  pioces.s,  l  or  example,  it  the  source  zone  concentration  "half- 
life"  is  10  years  and  the  initial  .source  zone  concentration  is  I  mg/L,  then  the 
source  zone  concentration  will  be  0.5  tng/1.  after  10  years,  and  0.2.8  mg/I, 
after  20  years. 

Note  lliat  the  assumption  that  di.ssolution  is  a  first-order  process  is  only  an 
approximation,  and  that  source  .ittenuation  is  best  described  by  first-order 
decay  when  concentrations  are  relatively  low  (<c  1  mg/I.).  For  more 
information  on  dissolution,  see  Newell  ct  al.,  (1944).  The  source  half-life 
IS  NOT  related  tii  lambda,  the  biodegradation  hall-lite  foe  dissolved 
constituents.  Lambda  is  used  to  calculate  the  .imouiit  of  biodegr,id.ition  of 
dissolved  organics  u/tcr  they  leave  the  .source  zone  and  travel  through  the 
plume  area  The  source  half-life  is  related  to  the  rate  of  dissolution 
ivcumng  111  the  source  zone,  anti  describes  the  change  in  .stiurre 
concentrations  over  time. 

•  The  BIOSCREEN  software  'utomatically  calculates  the  s<.)urce  zone 
concentration  hall-lite  it  the  user  enters  a  best  estimate  for  the  mass  ot 
liissolvable  organics  zone  (st'luble  organic  constituents  sorbed  m  the  soil, 
residual  NAI’Ls,  and  tree  product)  in  the  source.  The  half-life  ol  the 
di.ssolution  prcxess  c.in  b«'  approximated  it  one  knows  the  mass  ol 
dissolvable  organics  in  the  source /.one  (in  mg  or  kg),  the  How  rate  through 
the  source  zone,  and  the  average  concentration  ot  dissolved  organics  that 
leave  the  source  zone  The  equation  is  based  cn  integrating  the 
concentration  vs.  time  relationship  (tirst-order  decay)  and  using  the 
rel.itionship  that  the  ni,iss  in  the  sounv  zone  over  time  is  proportional  to 
the  source  cimceotration  over  time.  Ihis  yields  the  following  expression  tor 
the  halt-lite  ot  the  concentration  ot  dissolved  org.inics  in  the  source  zone 
(see  Appendi.x  A. ,5): 


t ,  „  -  (0.<)93  *  Mo  )  /  (Q  *  Co) 

lulf  wiunt  ''  ^  '' 


.  I 


MalHiff  nt  sDurci’ 
concentration  (yrs) 

Cin'uiulvvater  How  throu^la 
source  ^nne  (I,/yr) 

Effective  soun'e  /one  corv  . 
(obser\a’il  conf  entr.itioii  ♦  huule^; 
capacity  for  inst  ifact. 
assumption)  at  t  t)  (ir.^/I.) 

fVla.ss  ot  ilissf>lval)U‘  oT^anu's 
in  source /one  lit  f  0  (my,) 
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Description  (cont’d) 

Key  Questions; 

Why  are  there  two  source  half-lives  reported?  Note  that 
BIO:  CREEN  automatically  selects  the  correct  source  half-life  value 
depending  on  which  kinetic  model  is  being  used  (see  Which  Model 
Should  One  Use?  under  BIOSCREEN  Concepts). 

Two  .source  halt-lives  are  reported  by  the  model  in  the  .sciua'e  hall- 
life  cell:  the  smaller  num^  will  be  the  source  half-life  trom 
dissolution  if  Instantaneous  Reaction  kinetics  are  used,  and  the  larger 
value  will  be  for  No  Degradation  or  First-order  Decay  kinetics.  The 
first-order  decay  model  assumes  biodegradation  starts  immediately 
downgradient  of  the  source,  and  that  the  rate  of  dissolution  is 
reflected  bv  the  concentr  ition  of  dis.solved  organics  act  allv 
measured  in  monitoring  wells.  In  other  woras,  the  tirst-order  d  -cav 
model  assumes  Cq  is  equal  to  tlie  observed  source  concentration. 

The  instantaneous  reaction  model  assumes  biodegradation  is  occurrmg 
directlv  m  the  sourev  zone,  and  that  the  effective  source  zone 
concentration  C'o  is  equal  to  the  measured  concentration  in  the  stiurce 
zone  plus  any  "mi.ssing"  concentration  due  to  biodegradation.  For 
example,  it  the  source  zone  concentration  in  monitoring  wells  is  5 
nig/L,  and  the  biodegradation  capacity  is  10  mg/[.,  the  effective 
soune  concentration  Qi  ncentration  before  biodegraaation)  is  lb 

mg/I  .  In  other  word.s,  Co  is  equal  to  the  measured  .s<nme 
concentration  plus  the  biodegradation  capacity  provided  by  the 
electron  acceptor  concentration.  This  means  use  ot  the  instantaneous 
reaction  assumption  will  result  in  higher  dis.solution  rales  and 
shot  ter  source  lifetimes  (  see  Newell  et  «/.,  IW.S). 

Does  BIOSCRFEN  account  for  travel  time  away  from  the  declining 
source?  With  the  declining  source  option  in  BIOSCKEEN,  the 
concentration  tor  anv  location  and  any  time  is  calculated  using  a 
source  concentration  determined  by  the  first-order  decav  calculations 
shown  abiive  Hie  tinu*  used  to  determine  the  soun'c  concentration  is 
adjusted  to  ao  ount  lor  the  travel  time  between  the  soume  and 
measurement  point. 

I  n  "xample,  consider  the  casi'  where  a  declining  soune  term  is  ilscxI 
vvitli  a  source  lialt-lite  ot  1(1  years  and  a  solute  velocitv  ot  IIX)  It/vr. 

Ill  c.ilculati'  the  coiventralion  at  a  point  2(XX)  tt  aw.iy  at  time  bl) 

vears,  lilOSt  KlTiN  follow  -  these  stt'ps 

1)  (  .ilculc.tes  li.ive!  tune  trom  point  to  source:  2()(X)/  U)0  Jl)  ve.as 

2)  Subtracts  travel  lime  trom  simulation  time:  bO  yrs  -  2t)  vrs  -M  vrs 

3)  (  alcul.iies  source  decay  riH'lt.:  ksource -- l).(i'bb/(source  halt  liie) 

4)  C'alcuiates  source  cone,  at  1  -  iiiyr;  C.  „,  =  C„  exp' " 

Typical  Values 

1  to  1(),(X)()  years 

Source  of  Data 

C  alculateii  bv  model  from  soluble  mass  in  NAl'I..ind  soil  (see  hi'lovv), 
source  coiiLenlrations,  and  groundwater  velocitv 

How  to  Enter  Data 

C  .ill  iilated  liirectly  bv  model.  (  ba.nge  fiy  i  haiiging  soluble  mass. 

arametet 

*w  *»r 


Typical  Values 


Source  of  Data 


The  best  estimate  of  dissolvable  organics  in  the  source  zone  is 
obtained  by  adding  the  mass  of  dissolvable  organics  cn  soils,  free- 
phase  NAPf.s,  and  residual  NAPLs.  This  quantity  is  used  to 
estimate  the  rate  that  the  source  zone  concentration  declines.  Note 
that  this  is  an  experimental  and  unverified  model  that  should  be 
applied  with  care  (the  model  probably  underpredicts  removal  rate). 

F'or  ga.solme  or  JP-4  spills,  BTEX  is  usually  assumed  to  comprise  the 
bulk  ot  dis.solvable  organics  in  the  source  zone.  To  simulate  a 
declining  source,  use  the  method  described  below.  For  constant-source 
simulations,  either  enter  a  very  large  niimbtT  for  soluble  mass  in  the 
source  zone  (e.g.,  l,tXX).00{)  kg)  or  type  "Infinite". 


U,1  to  100,000  kg 


This  information  wilt  most  likely  come  from  cither: 

1)  Estimates  ot  the  ma.ss  of  spilled  fuel  (remember  to  convert  the 
total  mass  ot  spilled  fuel  to  the  dissolvable  mass;  for  example  HTliX 
represents  only  5-15%  ot  the  total  mass  ot  gasoliiu'). 

2)  Integration  of  maps  showing  contaminated  soil  zones  (data  in 
mg/kg)  .md/or  NAPl.  zonis  (usually  product  thickness).  The  user 
should  estimate  the  volume  of  contaminated  soil,  convert  to  kg  of 
contaminated  soil,  and  multiply  by  the  average  soil  concentration. 
To  make  the  estimate  more  accurate,  the  user  might  have  to  divide 
the  soil  into  different  zones  ot  soil  concentrations,  into  unsaturated 
\'s.  saturated  soil,  and/or  into  ditterent  depths.  (One  st.ind.ird 
approach  is  to  divide  into  a  vertically  averaged  iinsatui-.iled  /one 
map  and  a  vertically  averaged  saturated  zone  map.)  II  the  ust'r  is 
making  estimates  trom  NAPl  data,  remember  the  thickness  ot 
product  in  a  aquiter  IS  only  Ul-5l)'’i' of  the  actual  product  thickness  in 
the  well  (Hedient  ct  ul.,  l'>04) 

Note  tliat  the  data  is  to  b*‘  entereii  in  kg,  and  the  modi'l  will  convert 
the  results  to  estimate  the  source  halt-lite  An  example  is  pro\  nii'd 
below  assuming  a  bulk  density  ol  1  7  kg/I.  (e  g.,  l(X)  ft"  x  2(1  tl  '  2M.  ( 
l./tl' X  1  7  kg/I  \  mHi  mg/Kg  X  hi'"  kg/mg  ■  kg) 


SOLUBLE 

MASS 


100  iq.  ft  Drpth  20  ft 

Average  Soil  Cuticentration 
-  600  mg/Kg  BTEX 


Soil  Zone  2;  220  iq.  ft 


Depth  20  ft 


Average  Soil  Concentration 

-  SO  mg/Kg  BTEX 

'  Soil  Zone  };  400  iq.  ft  Depth  20  ft 

Average  Soil  Concentration 

-  10  mg/Kg  BTEX 


TOTAL  SOLUBLE  MASS 
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7.  FIELD  DATA  FOR  COMPARISON 


Units 

mg/1. 

Description 

These  parameters  are  eoni'entnitions  ot  dissolv'od  org.inies  in  wells 
near  the  s'enterline  of  the  plume.  These  data  are  usTtl  to  help 
calibrate  the  model  .md  are  displayed  with  model  results  in  the 
"Run  C'enterline'  option. 

Typical  Values 

Il.tXll  to  50  mg/T 

Source  of  Data 

Monitoring  wells  located  near  the  centerline  ot  tlu'  plume. 

How  to  Enter  Data 

1. liter  as  mans  or  .is  tew  ot  these  points  .is  needed.  1  he  data  are  usted 
only  to  hiTp  c.ilihr.ite  the  model  sx'hen  comparing  the  results  troni 
the'  centerline  optum.  Note  that  the  distance  from  source  values 
c.innot  be  changed:  use  the  closest  s  .due  possible 
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ANALYZING  BIOSCREEN  OUTPUT 

rhf  output '.Ikh'.’s  thr  Li-nturlinc  (lor  all  throe  kiiielu'  i  uhIoIs  at  the  saiiio 

tiiuo)  or  a^  an  array  (one  kmotio  model  at  a  time)  Note  that  the  results  are  .dl  tor  the  tune 
eutiTi'd  in  the  "Sinuil.ition  Time"  ho\^ 

Centerline  Output 

r  eiiterline  output  is  displayed  when  the  "Run  C  enterline"  Inilton  is  [Hesseil  ui  the  input  serei-n. 
Hie  i  i'iiterline  cnitput  seris'ii  shows  the  ar'era^e  eoneeniration  at  the  lop  o(  the  saturated  /one 
1/  III  alon^  the  c  enterline  ot  the  plume  (V  -Oi  C 'lu  kiiiy;  m ''Animate"  divides  the  simulation 
iiiti'  111  se[i, irate  tune  periods  .uul  shows  tlw  movement  i>l  tlie  plume  Tased  cii  the  three 
nil '‘■T  Kl'.TN  models  (red  ivc  degradation,  hlue-  tirsl-order  decsu',  green:  instant. ineous 
re.ietionl.  Note  th.lt  ,ill  c dneentr.itioiis  .ire  displ.ived  in  units  ol  mg,'  I  . 

Array  Output 

1  he  .irr.iv  output  1-  displ.n  isl  when  the  "Run  ,\rr,ie  "  inilton  is  pres.-ed  m  the  Input  sc ns'ii.  1  he 
User  is  .isked  to  selec  t  one  ol  the  lluee  model  types  (nc'  deglsicl.illon,  tus|  urdei  dec  .ly,  01 
instil  1 1. meous  re.ietionl  '  !>  er.ijihic  shows  ri’siills  m  .1  10  point  lonp,  he  point  wide  grid 
111  .iltei  the  iiuulelecl  .ire.i,  .ul|ti.t  the  Model  .Are.i  I  eilglli  .Hid  Width  p.ir.lllielels  m  the  input 
sc  i  i’en 

lo  SIC  the  plume  .irrav  th.lt  esceeds  .1  u  rt.im  l.irget  h-Vel  (sueh  .e.  .in  Ml  I.  or  risk-b.ised 
.  le.iiiup  Icwel),  rniei  the  i.irgi’t  le\'el  111  liie  i'ox  .iiiii  push  ■  Tint  1  Alt.i  ■  l.irget  I  >nlv  sei  lions  ol 
the  plume  exceeding  the  t.iiget  level  will  he  displ.ived  losee  .ill  the  d.ii.i  .igam,  push  "Tlol 
All  It.it.i"  Note  th.lt  lilUSi  Kl'l  N  . lilt. miatic  .illy  resets  this  hiittoii  to  "Plot  All  U.it.i''  when 
the  "Run  .Arms  ''  Inilton  1,  pressed  m  the  input  sereeu  An  .ipproxim.ite  ni.i  .s  h.il.inic  o. 
presented  on  the  .irr.iv  output  sc  reeii  .Is  liesi  riheci  helosx 
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Calculating  the  Mass  Balance 


The  model  calculates  the  total  amount  of  dissolved  contaminant  that  has  left  the  source  zone. 
If  the  source  is  an  infinite  source,  then  the  calculation  is  based  on  the  discharge  of 
groundwater  through  the  source  zone  (Darcy  velocity  for  groundwater  times  the  total  sciurce 
width  times  the  source  depth)  tunes  the  average  concentration  of  the  source  zone  (a  weighted 
average  of  concentrahon  and  source  length  for  each  of  the  different  source  zones)  times  the 
simulation  time. 

If  the  .soume  is  a  declining  source,  an  exponential  source  decay  term  is  ust.cl  to  estimate  the 
ma.ss  of  organics  that  have  left  the  source  zone  (see  Source  Data:  Varying  Concentrations 
Dver  Time).  Note  that  the  .source  decay  term  is  for  dissolution  of  soluble  organics  from  the 
source  zone  and  is  not  related  to  the  first-order  decay  term  for  the  dissolved  constituents. 

Note  that  the  total  mass  ui  the  plume  is  the  same  for  the  No  Degradation  and  First-order 
Decay  models  but  is  different  for  the  Instantaneous  Reaction  model.  The  source  zone 
dissolution  rate  is  calculated  to  be  much  higher  if  the  instantaneous  reaction  model  is 
selected.  The  Instantaneous  reaction  a.ssumes  that  active  biodegradation  reactions  ixrLiir  in 
the  source  zone,  and  that  the  observed  concentrations  of  organics  in  .stiurce  zone  monitoring 
wells  reflect  conditions  alter  biodegradation.  In  this  case,  the  actual  concentration  of 
organics  coming  oft  the  source  zone  is  Hjual  to  the  measured  concentration  plu.s  the 
biodegradation  capacity  of  the  upgradient  groundwater.  The  resulting  higher  effective 
dissolution  rate  equ.ites  to  a  greater  amount  of  mass  leaving  the  source  area,  leading  to 
ditterent  mass  values  tor  the  Instantaneous  Reaction  model. 


ActuiU  Plume  Masi(kg)-jjfr! 

HIObCKld'.N  cali  ulati's  the  m<is:.  ot  organics  in  the  5x10  p/lume  array  for  the  three  models: 

1 )  No  Degradation  2)  1st  Order  [.)ecay  3)  Instantaneous  Reaction 

The  mass  is  calculated  by  assuming  that  each  point  ropresents  a  cel  -qual  to  the  incremental 
width  .iiul  length  (except  hir  the  tir.f  column  which  is  a.ssiimixi  to  bx'  half  as  long  as  the 
otfuT  columns  bee, ui.se  the  source  is  .issunuxi  to  'x'  in  the  miikile  of  the  cell).  The  volume  ot 
.ittected  groundw.iter  in  e.ich  cell  is  calculated  hv  multiplying  the  ,irea  ot  e.ich  cell  bv  the 
s<iuri'e  depth  and  by  porosity  (the  inasi^  b.ilance  calculation  assumes  2-D  transport).  1  he  mass 
ot  org.iiiii  s  in  each  cell  is  then  determined  hv  nuiltiplving  the  volume  of  groundwater  hv  the 
conccnti.ition  and  tlicn  bv  the  ret,ird,Uion  tactur  (to  account  tor  sorlu'd  constitiiciits). 

How  BIOSCREEN  Estimates  Actual  Plume  Mass  for  BiodegraH ition  Models 

It  the  mass  ot  orgiiiiii  s  in  the  5x10  plume  ,irr,\y  is  witfiin  50";,  to  i  >0'!;,  ot  the  mass  ot  org.inics 
th.it  have  left  the  source  (see  box  .ihove),  then  two  values  are  e.ileiil.ited: 

"li  Biodegraded,  1st  order  decay  '  (Plume  Mass,  1st  oriler  ilecay)  *  IIX)  /  (Plume  mass,  ixi 
huuieg) 

Huidegraded,  inst.  re.u't.  (Plunu'  Mass,  iiist.  react)  *  100  /  (Plume  m.iss,  no  biodeg) 

these  pereent.iges  are  inultifilied  .ig.iinst  the  Plume  Mass  if  No  Biodegradation  Value  (first 
box)  til  estim.ite  the  aetii.il  plume  mass  tor  the  tvvo  hioi.Iegr,ui.ition  models.  If  the  No 
I  tegr.id.ition  model  lias  been  selected,  there  is  no  bloiiegnul.ituiii,  ,ind  the  Actual  Plume  Mass 
(second  !iox)  will  ei|iMl  ihe  Plume  M.is.s  if  No  Bio’degrad.ition  (first  box). 
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If  BIOSCREEN  Says  “Car.’t  Calc" 

If  the  mass  of  organics  in  the  plume  does  not  fall  within  50%  to  150%  of  the  mass  of  organics 
that  have  left  the  source  (first  box),  then  the  model  concludes  that  the  modeled  area  (see 
Input  Screen,  Section  5:  General  Data)  is  not  sized  correctly  to  capture  enough  mass  in  the 
5XiO  array  and  writes  "Can't  Calc"  in  the  box.  The  user  is  en.  ouraged  to  adjust  the  modeled 
length  and  width  to  capture  most  of  the  No  Degradation  plume  in  the  5x10  array.  In 
addition,  sometimes  source  conditions  with  variable  concentrahons  and  widths  (see  input 
screens)  can  make  it  difficult  to  accurately  capture  the  plume  mass.  If  the  user  has  problems 
obtaming  a  mass  balance  even  after  changing  the  modeled  area,  change  the  source  term  to  a 
single  source  zone  (instead  of  3  or  5  zones)  to  improve  the  accuracy  of  the  mass  balance. 

If  problems  still  exist,  ensure  that  the  vertical  dispersivity  term  (Section  2  on  the  Input 
Screen',  is  set  to  0  (the  default  value).  The  mass  balance  calculations  are  less  accurate  for 
three-dimensional  simulations. 


.Plume '.Mass  Remoyed  by^Biode2ra^tioi^(lcj^‘^^^ 

mmm 

An  estimate  of  the  rriass  of  conto'  .nanls  that  are  biodegraded  is  provided  in  BIOSCREEN. 
The  model  subtracts  the  Actual  Plume  Mass  (second  box)  from  the  Plume  Mass  if  No 
Biodegradation  (first  box).  For  the  No  Degradation  model,  the  first  box  equals  the  second 
box,  and  Plume  Mass  Removed  by  Uiodeg  is  zero.  For  the  other  two  cases,  the  2  boxes  will 
differ,  and  the  amount  of  biodegradation  will  be  calculated.  The  value  beneath  the  third 
box  shows  the  %  of  organics  tha*-  have  left  the  source  and  have  been  biodegraded. 


Change;  in  Elec^n,  Acceptor/Bypraduct^Masm; 


nOSCREEN  u.scs  the  Plume  .Mass  Removed  by  Biodegradation  to  back-calcuUte  the  amount 
of  measurable  electron  acceptors  consumed  and  the  amour.t  of  mea.surable  metabolic  by¬ 
products  that  have  been  produced. 

For  example,  the  amount  of  oxygen  consumed  is  calculated  by: 

Oxygen  Consumed  (kg)  =  (Plume  Mass  Removed  by  Biodeg)  ’  (Delta  Q2/Util.  Fad.) 

(  Biodeg.  Capacity) 

(see  Biodegradation  Capacity  section  to  see  how  this  term  i.s  calculated) 

Note  that  the  total  .sum  of  consumed  electron  acceptors  does  not  equal  the  Plume  Ma.ss 
Removed  by  Biodegradation.  This  is  because  the  stoichiometry  of  the  biodegradation 
reactions  do  not  represent  a  1:1  relationship  betwetm  the  ma.ss  of  hydrocarbon  and  electron 
acceptor  consumed  (see  Utilization  Factor  sec'ion) 


Oripi  al  iMuS;  in'Sourej^ 


Equal  to  the  Soluble  Ma.ss  in  NAPE  and  Soil  entered  by  the  user  cn  the  Input  Screen.  If  the 
user  has  selected  an  "Infinite"  mass  to  simulate  a  non-declining  source,  this  box  will  show 
"Infinite." 


The  amount  of  mass  n  maiiung  in  the  .source  zone  at  the  end  of  the  simulation  pei  iod  is 
calculated  and  displayed  in  this  box.  This  calculation  is  performed  as  follows: 

(Mass  in  the  .Source  Now)  - 

(Original  Mass  in  Source)  (Actual  Plume  Ma.ss  e  Plume  Ma.ss  Removed  by  Biodeg) 
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If  the  mass  of  orgzinics  in  the  plume  falls  within  50%  to  150%  of  the  mass  of  organics  that 
have  left  the  source  (fir'  i  box),  then  the  model  concludes  the  modeled  area  (see  Input  Screen, 
Section  5;  General  Data)  is  appropriately  sized  to  estunate  the  volume  of  the  plume.  In  this 
case  BIOSCREEN  counts  the  number  of  cells  in  the  5  x  10  array  with  concentration  values 
greater  than  0,  and  mnlhplies  this  by  the  volume  of  groundwater  in  each  cell  (length  *  width 
*  source  thickness  ’  porosity). 

If  the  user  wishes  to  estimate  the  volume  of  the  plume  above  a  certain  target  level,  enter  the 
target  level  in  the  appropriate  box  and  pre.ss  the  appropriate  mo  iel  to  display  the  result 
(No  Degradation,  1st  Order  Decay,  or  Instantaneous  Reaction). 

Note  that  the  model  does  not  account  for  the  effects  of  any  vortical  dispt'i'sion. 


Flowrate  of  .Water  Through 

Using  the  Darcy  velocity,  the  source  thickness,  and  the  source  width,  BIOSCREEN  calculates 
the  rale  that  clean  groundwater  moves  through  the  .source  /one  where  it  will  pick  up 
dissolved  hydrocarbons.  Note  that  the  groundwater  Darcy  velocity  is  equal  to  the 
groundwater  seepage  velocity  multiplied  by  porosity. 
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BIOSCREEN  TROUBLESHOOTING  TIPS 
Minimum  System  Requirements 

The  BIOSCREEN  model  requires  a  computer  system  capable  of  running  Microsoft*  Excel  S.O  for 
Windows.  Because  of  the  volume  of  calculations  required  .o  process  the  numerical  data 
generated  by  the  model,  GSI  recommends  ruruiing  the  model  on  a  system  equipped  with  a  486  DX 
or  higher  processor  running  at  66  MHz  or  faster.  A  minimuiTi  of  8  Megabytes  of  systciu  memory 
(RAM)  is  strongly  recommended. 

The  model's  input  and  output  screens  are  optimized  for  display  at  a  monitor  resolution  of 
640x480  (Standard  VGA).  If  you  are  using  a  higher  resolution,  for  example  800x600  or  1024x768, 
see  Changing  tiie  Model's  Display. 

For  best  results,  Start  Excel  and  Load  the  BSCREEN.XLS  file  from  the  File  /  Open  menu. 


Spreadsheet-Related  Problems 

The  buttons  won't  work:  BIOSCREEN  is  built  in  the  Excel  spreadsheet  enviromnent,  and  to 
enter  data  one  must  click  anywhere  outside  the  cell  where  you  just  entered  data.  It  you  can  see 
the  numbers  you  just  entered  in  the  data  entry  part  of  Excel  above  the  spreadsheet,  the  data  has 
not  yet  been  entered.  Click  on  anotlier  cell  to  enter  the  data. 

####  is  displayed  in  a  number  box:  The  cell  format  is  not  compatible  with  the  value,  (o.g.  the 
number  i.s  too  big  to  fit  mto  tlie  window).  To  fix  this,  select  the  cell,  pull  down  the  format  menu, 
select  'Cells"  and  click  on  the  "Number"  tab.  Change  the  format  of  the  cell  until  the  value  is 
visible.  If  the  values  still  cannot  be  read,  select  tlie  format  menu,  select  "Cells"  and  dick  on  the 
"Font"  tab.  Reduce  the  font  size  until  the  value  can  be  read. 

#DIV/0!  i.s  displayed  in  a  number  box:  The  most  aimman  cause  of  this  problem  is  that  some 
input  data  ,ire  missmg.  !n  some  cases,  entering  a  zero  in  a  box  will  cause  this  problem.  Double 
chock  111  make  certain  that  all  of  tlie  input  cells  required  tor  your  run  have  data.  Note  that  for 
vertical  dispersivity,  BIOSCKEEN  will  convert  a  "d"  into  the  data  entry  cell  into  a  very  low 
numiier  (1x10  ■")  to  avoid  #DIV/0!  errors. 

There  once  were  formulas  in  some  of  the  boxes  on  the  input  screen,  but  they  were  dccidentally 
overwritten:  Click  on  the  "Restore  Formulas  for  Vs,  Dispersivities,  K,  and  lambda"  button  cn 
the  bottom  right-hand  side  of  the  input  screen.  Note  that  this  button  will  also  restore  the 
formulas  that  make  the  Source  Width  and  Source  Concentrations  for  source  zones  4  and  5  equal 
to  source  zones  2  and  1,  respectively. 

The  graphs  .spem  tn  nii  .ye  around  and  change  size:  Thi;.  is  a  feature  of  Excel.  When  graph  scales 
are  altered  to  accommodate  different  plotted  data,  the  phy.sical  size  of  the  graph',  will  change 
slightly,  sometimes  resulting  in  a  graph  that  spreads  out  over  the  fixed  axis  legends.  You  can 
manually  resize  the  graph  to  make  it  look  nice  again  bv  double-clicking  in  the  graph  and 
resizing  it  (refer  to  the  Excel  U.ser's  Manual). 

Comrtio.n  Error  Messages 

Unable  to  Load  Help  File:  The  luo:.*  nimniin  ermr  message  encounlured  with  BRXSCRFHN  is 
the  message  "Unable  to  Dyreii  Help  I'ile"  after  ciu 'r-iiig  m  a  Help  button  Depeiidmg  m  the 
version  of  Windows  vou  are  using,  you  may  get  an  I'm  el  I  )i,iiiig  Box,  -i  Wiiiilows  I  .'ialog  lio.x,  or 


you  may  see  Windows  Help  load  and  display  the  error.  This  problem  Ls  rela  'd  to  the  ease 
with  which  the  Windows  Help  Engine  can  f^d  the  datafile,  BIOSCRN.HLP.  Here  are  some 
suggestions  (in  decreasing  order  of  preference)  for  helping  WinHelp  find  it: 

•  If  you  are  fortunate  enough  to  be  asked  to  find  the  requested  datafile,  do  so.  It's 
called  BIOSCRN.HLP,  and  it  was  installed  in  the  same  directory /folder  as  the 
BIOSCRN.XLS  file. 


Use  the  File/Open  menus  irom  within  Excel  instead  of  double-clicking  on  the 
filename  or  Program  Manager  icon  to  open  the  BIOSCRN.XLS  file.  This  sets  the 
"current  directory"  to  the  directory  containing  the  Excel  file  you  just  opened. 

Change  the  Winlleip  rail  in  the  VB  Module  to  "hard  code"  the  directory 
information.  That  wa}  ,  the  file  name  and  its  fuH  path  will  be  explicitly  pa.ssed  to 
WinHelp.  Hints  for  doing  this  are  in  the  VBA  module.  Select  the  Macro  Mod.'ili-  tab 
and  search  for  the  text  "Helpfilo". 


•  Asa  last  resort,  you  can  add  the  BIOSCREEN  directory  to  your  path  (located  in  your 
AUTOEXEC.BAT  file),  and  this  problem  wiU  be  aired.  You  will  have  to  reboot  your 
machine,  however,  to  make  this  work 

The  BIOSCREEN  system  was  designed  to  be  used  m  a  PC  with  Windows  configured  to  a 
standard  VGA  resolution  of  640x480  pixels.  If  you  are  using  a  larger  monitor  and  your  video 
resolution  is  set  to  800x600  pixels  or  greater,  you  will  need  to  change  the  zoom  factor  in  the 
Visual  Basic  code. 

In  the  first  three  Lines  in  the  Macro  Module  of  the  BIOSCREEN  spreadsheet,  change  th" 
number  after  the  equals  sign  in  the  following  line: 

Const  ZoomValue  =  65 

If  your  display  resolution  is  standard  VGA  (640x480),  use  65  tor  the  zoom  value.  If  your 
resolution  us  800x600,  use  a  zoom  value  of  82.  If  your  re.solubon  is  not  640x480  or  800x600,  if  your 
video  pertormance  is  senously  degraded,  or  if  you  experience  display  problems,  you  may  need  I  > 
change  your  video  resolution  (see  the  online  help  for  Windows  Setup  or  consult  your  Windows 
installation  manuals)  and  experiment  dth  other  values  for  ZoumValue. 
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APPENDIX  A.  I  DOMENICO  ANALYTICAL  MODEL 


The  C)omenico  (1987)  analytical  model,  used  by  BIOSCREEN,  is  designed  for  the 
multidimensional  transport  of  a  decaying  coritaminant  species.  The  model  equation,  boundary 
conditions,  assumptions,  and  limitations  are  discussed  below. 

Domenico  Model  with  Inrtantaneoui  Reaction  Superposition  Algorithm  ’  '' W,  .  ' 

.■Sv  ^  .  C(x.y.o.f)  _  1  _ J  XI, 


(Oi+BC)  8  [a,2\  ^  ‘  ) 

|x  -  vt(l  +  /  v)'”  j 


f  s\7T  i  vW2 

[2(a,x)  J  2(a,x) 


(Z)  (-71 


BC  =zT 


Definitions _ _ 

BC  DiodcgradaUon  capacity  (m^/I ) 

C(x,y.z,t)  Conccixtration  at  di*»Unce  x  dowixstream  of 

source  and  dL^Unre  y  off  crtxterlme  of  plume  at 
tune  t  (in^/l.) 

C.'^  Conccntratum  in  ^Kiurce  /tone  (mn/l.) 

Cn  C  oncentrdtinn  u\  Source  Zone  atl'-O  (mg/L) 

*  Distance  dov/ngradinit  of  source  (ft) 

y  DisUiice  from  plume  centerline  <if  v)«rce  (ft) 

J  Instance  from  surface  to  measurement  point 

(assumed  to  be  0;  conemtrarjon  l»  always 
rtViumed  to  be  at  top  of  water  table) 

t'(iM),^  ConcL»iiiralJon  of  electron  acceptor  n  in 
groundwiiler  (itik/I.) 


Ul'n  Utilization  factor  for  electron  acceptor  n  (i.e-,  mass  ratio 
of  electron  acceptor  to  hydrocarbon  consumed  m 
biodei^adabon  reaction) 

Loiiptudmal  groundwater  dispersivily  (ft) 

Uy  Transverse  groundwater  dLspersivity  (ft) 
u,  VerticaJ  groundwater  dispersmty  (ft) 

0^  Effective  St>il  Furo.sity 

X  First-Order  Degradation  Rate  (day'M 

u  Groundwater  Seepage  Velrvity  (ft/yr) 

K  Hydraulic  Conductivity  (ft/yr) 

R  Constituent  retardation  factor 

i  Hydraulic  (.jradienl  (cm/ cm) 

Y  Source  Width  (ft) 

Z  Source  I,)epth  (ft) 


Tho  initial  conditions  arc: 

1)  c(x,  y,  z,  0)  -  0  (Initial  concentration  =  0  for  x,  y,  z,  >  0) 

2)  c(l),  Y,  Z,  ())  =  C„  (Source  concc-iiiiation  tor  each  vertical  plane  source  =  Cg  at  time  0) 


Tluj  1‘ey  assuuiptioiis  Ul  the  utodel  arci 

1)  The  aquifer  and  flow  field  are  homogenonous  and  isotropic. 

2)  The  groundwater  velocity  is  fast  enough  that  molecular  m  the  dispersion 

terms  can  be  ignored  (u\ay  not  be  arinr»>-;;^n;  lor  simulation  of  transport  through 
clays). 

'■)  /lusorption  is  a  reversible  process  represented  by  a  linear  isothenn. 
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The  key  limitaticas  to  the  model  arc; 


1)  The  model  should  rot  be  applied  where  pumping  systems  create  a  complicated  flow 
field. 

2)  The  model  should  not  be  applied  where  vertical  flow  gradients  affect  conf^nuiiant 
transport. 

3)  The  model  should  not  be  applied  where  hydrogeologic  conditions  change 
dramatically  over  ‘die  simulation  domain 


The  most  important  modifications  to  the  original  Domenico  model  are; 


0 


1)  The  addition  of  "layer  cake"  source  terms  where  three  Domenico  models  are 
superimposed  one  cn  top  of  another  to  yield  the  .S-  source  term  used  in  BIOSCKEEN 

■  (see  Connor  i  t  al.,  1994;  and  the  Source  Width  description  in  the  BIOSCREEN  Data  9 

Entry  Section). 

2)  Addition  of  the  instantaneous  reaction  term  using  the  superposition  algorithm  (see 
Append!);  A.2,  below).  For  the  instantaneous  reaction  assumption,  the  source 
concentration  is  .issumed  to  be  an  "effective  source  concentration"  (Coe)  equal  to  the 
observed  concentration  ui  the  source  zone  plus  the  biodegradation  capacity  (see 

"Source  Concentration"  on  the  BIOSCREEN  Data  Entry  section).  ® 


i 


» 
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APPENDIX  A  7  'SSTANrANEOUS  REACTION  -  SUPERPOSITION  ALGORITHM 

Early  biodegradatiorr  research  focused  on  the  role  of  dissolved  oxygen  in  controbing  the  rate  of 
biodegradation  in  the  subsurfaci'  (Borden  et  al.,  1986;  Lee  et  al,  1987).  Because  microbial 
biodegradation  kinetics  are  relatively  fast  in  comparisoxr  to  the  rate  of  oxygen  transport  in  the 
groundwater  flow  system,  Borden  demonstrated  that  the  biodegradation  process  can  be 
simulated  as  an  instantaneous  reaction  between  the  organic  contaminant  and  oxygen.  This 
simplifying  assumption  was  incorporated  into  the  BIOPLUME  I  numerical  model  which 
calculated  organic  mass  loss  by  supciposition  of  background  oxygen  concentrations  onto  the 
organic  contaminant  plume.  In  BIOPLUME  II,  a  dual-particle  mover  procei'ure  was 
incorporated  to  more  accurately  simulate  the  separate  transport  of  oxygen  and  organic 
contaminants  within  the  subsurface  (Rifai  et  al,  1987;  Rifai,  et  al,  1988). 

In  most  analytical  modeling  applications,  contaminant  biodegradation  is  estimated  using  a 
first-order  decay  equation  with  the  biodecay  ha'f-life  values  determined  from  research 
literature’  or  site  data.  However,  by  ignoring  oxygen  limitation  effects  such  first-order 
expressions  can  significantly  overestimate  the  rate  and  degree  of  biodegra''  ition,  particularly 
within  low-flow  regimes  where  the  rate  of  oxygen  exchange  in  a  groundwater  plume  is  very 
slow  (Rifai,  1994).  As  a  more  accurate  method  of  analysis,  Newell  recommended  mcorporation 
of  the  concept  of  oxygen  superposition  into  an  analytical  model  (Connor  et  al.,  1994)  in  a  manner 
similar  to  that  employed  in  the  original  BIOPLUME  model  (Borden  et  al.  1986).  By  this 
method,  contaminant  mass  concentrations  at  any  location  and  time  within  the  flow  field  are 
corrected  by  subtracting  1  mg/  L  organic  mass  for  eacli  3  mg/L  of  background  oxygen,  in  accordance 
with  the  instantaneous  reaction  assumption.  Borden  et  al  (1986)  concluded  this  simple 
superposition  technique  was  an  exact  replacement  for  more  sophisticated  oxygen-limited 
models,  as  long  as  the  oxygen  and  the  hydrocarbon  had  the  same  transport  rates  (e.g., 
retardation  factor,  R  =  1). 

Ill  their  original  work,  Borden  et  al.  (1986)  noted  that  for  highly  sorptive  rontaminaius  tne 
oxygen-superposition  method  might  erroneously  characterize  biodr j,iauuoon  due  to  the 
[liffering  transport  rates  of  dissolved  oxygen  and  the  Cigaiuc  contaminant  within  the  aquifer 
matrix.  However,  as  demonstrated  hy  Cuiuioret  at.  (1994),  the  oxygen  superposition  method 
and  BICPLLIME  If  td  j.d  particle  transport)  are  in  reasonable  agreement  for  contaminant 
retan-i  tactors  as  high  as  6.  Therefore,  the  superposition  method  can  be  employed  as  a 
reasonable  approximation  in  BIOSCKEEN  regardless  of  contaminant  sorption  characteristics. 

lilOSCREHN  employs  the  same  supeiposition  approach  fi'r  all  ot  the  aerobic  and  anaerobic 
biodegradation  reactions  (based  on  evaluafio.i  of  Oj,  NO,,  SO,,  Eo^',  and  CH,).  Based  cn  work 
report'd  by  Newell  ct  al.  (1995),  the  anaerobic  reactions  (nitrate,  ferric  iron,  and  sulfate 
reduction  and  methanogenesis)  are  amenable  to  simulation  using  the  instantaneous  reaction 
assumption.  The  general  approach  is  pre.sented  below: 
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Based  on  the  biodegradation  capacity  of  electron  acceptors  present  in  the  groundwater  system 
this  algorithm  will  correct  the  non-decayed  groundwater  plume  concentrations  pred-cUj  by  the 
Domenico  model  (Appendix  A.l)  for  the  effects  of  organic  constituent  biodegradation. 

4  To  summarize: 

„.iginal  BIOPLUME  model  (Borden  efal.  1986)  used  a  superposition  method  to  simulate 
the  fast  or  "instantaneous"  reaction  of  dissolved  hydrocarbons  with  dissolved  oxygen  in 
groundwater. 

Borden  et  at.  (1986)  reported  th.i.t  tliis  version  of  BIOPLUME  was  mathematically  exact  for 
the  case  where  the  retardation  factor  of  the  contiiminant  was  1.0. 

Rifai  and  Bedient  (1990)  developed  the  BIOPLUME  II  model  with  a  dual-particle  tracking 
routine  that  expanded  the  original  BIOPLUME  model  to  handle  contaminants  with 
retardation  factors  other  than  1.0,  in  addition  to  otlier  improvements. 

Connor  el  al.  (1994)  compared  the  superposition  method  with  the  more  sophisticated 
BIOPLUME  II  model  and  determined  that  the  two  approache.s  yielded  very  similar  results 
for  readily  biodegradable  contaminants  with  retardation  factors  between  1.0  and  6.0. 

BIOSCREEN  was  developed  using  the  supierposition  approach  to  simulate  the 
"instantaneous"  reaction  of  aerobic  and  .inaerobic  reactions  in  groundwater.  The  biodegra¬ 
dation  term  in  BIOSCREEN  is  mathematically  identical  to  the  approach  used  in  the 
original  BIOPLUME  model.  This  mathem.atical  approach  (supicrposition)  matches  the  more 
sophisticated  BIOPLUME  II  model  very  closely  for  readily  biodegradable  contaminant 
retardation  factors  of  to  6.0.  DIOSCREEN  simulations  using  file  instantaneous  reaction 
assumption  at  sites  with  retardation  factors  greater  than  6.0  hould  performed  with 
caution  and  vf*rifipH  rising  a  more  sophisticated  model  such  as  BIC9PI.Uh4E  ill. 

« 


1) 

2) 

3) 

4) 


5) 


Iftlo'S!  '"'***’  ™  d«c«y 

(but  with  sourc*  ion« 
coflcantrMlon  (qual  to 
im.uurwd  loump  zorw 
coiKwntmlon  ' 
blaa.'qradalion  capacity  SC  ) 


i-TriuLT'l  Subtract  Blotiagradatlon 
l-:*y.£il  Capacl^(eC)fromNo 
Decay  Consantratlona 


I  PraCIct  blodagnKtad  pluina 

— - - 1  concaniratlons  aaauming 

Inatantanaoui  RaacUon 
Auuinptlon 


daily 


c.-c_ . 
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APPENDIX  A.3  DERIVATION  OF  SOURCE  HALF-LIFE 


Purpose:  Determine  the  source  half-life  relationship  used  in  BIOSCREEN  (see  Souit'e 

Half-Life  discussion  in  BIOSCREEN  Data  Entry  Section,  pg  30), 

Given:  1)  There  is  a  finite  amount  of  soluble  organic  CMiipounds  in  source  zone  (the 

area  with  contaminated  soils  and  either  free-phase  or  residual  NAPL. 

2)  These  organics  dissolve  slowly  fresh  groundwater  passes  through  source 
zone.  Assume  the  change  in  mas<;  due  to  dissolution  can  be  approximated  as 
a  first  order  proces*;: 

M(f)  =  Mo  ^ 

Procedure:  1)  Calculate  initial  mass  of  dissolvable  organics  m  source  zone,  M,, 

2)  Determine  initial  source  concentration  from  monitoring  well  data,  C,, 

3)  Apply  conservation  of  mass  to  a  control  surface  containing  source  zone. 

4)  Set  the  expressions  for  mass  at  time  t  ^  Q  based  on  dissolution  and 
conservation  of  mass  equal  to  each  other  and  solve  for  an  expression 
describing  the  concentrahon  at  time  t  ^0. 

5)  Apply  initial  conditions  for  concentration  at  time  t=0  and  solve  for  the  first 
order  decay  constant,  k,. 

Assumptions:  1)  Groundwater  flowrate  is  constant,  (2(f)=(3i) 

2)  Groundwater  flowing  through  the  source  zone  is  free  of  organic  compounds. 
This  implies  that  no  mass  is  added  to  the  system,  only  dissolution  occurs. 


(!) 


Calculations:  1/  Calculate  initial  mass  of  dissolved/soluble  organic  ciimpiund,  M„  by  using 
procedure  described  under  "Soluble  Mass  in  NAPL,  Soil"  page  in 
BIOSCREEN  Data  Input  section. 

2)  neterrnme  initial  concentration,  C„  of  organic  compxmnd  in  groundwater 
leaving  the  source  zone.  This  may  be  a  spatial  average,  maximum  value,  or 
other  value  representative  of  the  groundwater  concentration  leaving  the 
source  area  (Note  that  for  the  instantaneous  reaction  assumption,  C„equdl . 
the  concentration  observed  in  monitoring  wells  plus  tiie  biodegradation 
capacity  to  account  for  rapid  hioHegradation  reactions  in  the  .source  zone. 
See  "Soluble  Mass  in  NAPL,  Soil"  page  in  BIOSCREEN  Data  Input  section). 

ai---l))  =  Co 

.3)  Apply  conservaten  of  mass  to  a  control  surface  that  contains  the  source 
zone.  The  mas.s  present  in  the  source  zone  at  time  t  ^  0  is  the  initial  mass 
plus  the  change  in  mass 

M(t)  M„+I{.^It  Q(t)C(t)dtdA 


(2) 


(3) 


4S 
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DERIVATION  OF  SOURCE  HALF-LIFE,  Cont'd 


Applying  the  assumptions  equation  (3)  simplifies  to 


« 

M(t)^  Mq  -  fi  Qo  at)  dt 

(4) 

• 

4) 

Set  the  two  expressions  for  mass  of  organic  compound  in  the  source  zone  a  t 
time  t  2  0  (equations  (1)  and  (4))  equal  to  each  other  and  solve  for  an 
expression  describing  the  concentration  leaving  the  source  zone. 

M„  e<‘  =  Mo -It  Qo  at)  dt 

(5) 

(6) 

• 

QQat)  =  l(.  M„e<^ 

(7) 

C(r)  = - -e 

Qo 

(8) 

• 

5) 

Apply  the  initial  condition  for  concentration  leaving  the  siiurce  zone  at 
time  f=0,  eqn  (2)  to  the  expression  for  at),  cqn  (8)  aiid  solve  for  the  first 
order  decay  coefficient,  k, 

r  -  ^0 

• 

0, 

(9) 

j.  _  QoQ-o 

Mo 

(11) 

Summary:  The  decay  coefficient  for  the  source  zone  in  BIOSCKEFN  is: 


'  Mo 


The  expression  for  mass  at  any  time  t  ^0  is: 

M(t)= 

Similarly  the  expression  for  source  zone  concentration  any  time  t  2  0  is: 
C<t)=  C„  <• 


Acknoiuledjimenl^-  Urifftnal  derivation  developed  lip  C.  Newell.  Detailed  derivation  developed  by  Xtaiiiiiiiiy 
l.iu,  Anthony  Holder,  and  Thoma.'i  Reever. 


« 
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APPENDIX  A.4  DISPERSIVirr  ESTIMATES 


Dispersion  refers  to  the  process  v/hereby  a  plume  will  spread  out  in  a  longitudinal  direction 
(along  the  direction  of  groundv^ater  flow),  transversely  (perpendicular  to  groundwater  flow), 
and  vertically  downwards  due  to  mechanical  mixing  in  the  aquifer  and  chemical  diffusion. 
Selection  of  dispersivity  values  Ls  a  difficult  pnxress,  given  the  impracticability  of  measuring 
dispersion  in  the  field.  However,  dispersivity  data  from  over  50  sites  has  been  compiled  by 
Gelhar  et  al.  (1992)  (see  figures  A.l  and  A. 2,  next  page). 

The  empirical  data  indicates  that  longitudinal  dispersivity,  in  units  of  length,  is  related  to 
scale  (distance  between  source  and  measurement  point;  the  plume  length;  Lp  in  BIOSCREEN). 
Gelhar  cf  a!.  1992)  indicate  there  is  a  considerable  range  of  dispersivity  values  at  any  given 
.scale  (on  the  order  of  2  -  3  orders  of  magnitude),  2)  suggest  using  values  at  the  low  end  of  the 
range  of  possible  dispersivity  value.s,  and  .3)  caution  against  using  a  single  relationship  between 
scale  and  dispersivity  to  estimate  dispersivity.  However,  most  modeling  studies  do  start  with 
such  simple  relationships,  and  BiOSCREEN  is  programmed  with  some  commonly  used 
rel.itionships  representative  of  typical  and  low-end  dispersivities: 


•  Longitudinal  Dispersivity 
Alpha  X  ^^  2H  O  S.t 


•  Transverse  Dispersivity 

Alpha  y  -  0.10  alpha  x 

•  Vertical  Dispersivity 

Alph.i  /  -  very  low  (i.e.  1 


(  )ther  coinmonlv  used  relationships  include: 


Alph.i 

X 

-  0.1  l.p 

(I’ll kens  and  Crisak.  19S1) 

Alph.i 

V 

jlph.i  X 

(ASTM.  1995)  (i:i'A.  I9B(,) 

Alph.i 

/ 

O.OS  alpha  x 

(A.^TM,  1995) 

Alpha 

/ 

alplia  X  tu  0.1  alpha  x 

(KPA.  I98ti) 

[  Ef  3  (Xu  and  l.'cfcjtfiii,  1995) 

\  .1.28  J 

f.,  in  Jt) 

( llasfJ  on  fiiy/i  reliabiliti/  B  ^ 

IH'iiils  from  Ctlhar  rl  ill.,  1992) 

1  e-99  ft)  (Based  on  conservative  eslimati 


The  BICISGRHEN  input  .siTeen  includes  Excel  fonriulas  to  cstiinate  diapersivitie.s  fiOiii  .scale. 
BIOSCKEEN  u-ses  the  Xu  and  Eckstein  (1995)  algorithm  for  estimating  longitudinal 
dispersivities  because  1)  it  provides  lower  range  estimates  of  dispersivity,  especially  for  large 
values  of  Lp,  and  2)  it  was  d(;veloped  after  weighting  the  reliability  ol  the  various  ficKl  data 
comjiili'il  by  Gelhar  et  a!..  (1992)  (see  Figure  A.l).  BIOSCKEEN  also  employs  low-end 
estimates  for  transverse  and  vertical  ili.spersivity  estimate.--  ,'>.10  alpha  x  and  0,  respectively) 
hccause:  1)  these  relationships  better  fit  observed  field  dit-’  reported  by  Gelhar  ft  al.  to  have 
high  reliability  (see  Figure  A. 2),  2)  Gelhar  t’f  al.  recommend  u.se  of  values  in  the  lower  range  of 
the  ofiserved  d.it.i,  and  3)  better  results  were  realized  when  calibrating  BIOSCREEN  to  actual 
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field  sites  using  lower  dispersivities.  The  user  can  override  these  formulas  by  directly  entering 
dispersivity  values  in  the  input  sceen  cell. 

Note  that  the  Domenico  model  and  BIOSCREEN  are  not  formulated  to  simulate  the  effects  of 
chemical  diffusion.  Therefore,  contaminant  transport  through  very  slow  hydrogeologic  regimes 
(e.g.,  clays  and  slurry  walls)  should  probably  not  be  modeled  using  BIOSCREEN  unless  the 
effects  of  chemical  diffusion  are  proven  to  hie  insignificant.  Domenico  and  Schwartz  (1990) 
indicate  that  chemical  diffusion  is  small  for  Pcclct  numbers  (seepage  velocity  times  median 
pore  size  divided  by  the  bulk  diffusion  coefficient)  greater  than  100. 


Figure  A.I.  Longiiudinal  dispcrsivity  vs.  scale  data  reported  by  Gcihar  ct  ol  (1992). 
Dau  includes  Gelhar's  re?nalysls  of  several  dispersivity  studies.  Size  of  circle 
represents  general  reliability  of  dispersivity  estimates.  Location  of  10%  of  scale 
linear  relationship  plotted  as  dashed  line  (Pickens  and  Grisak.  1981).  Xu  and 
Eckstein’s  regression  (used  in  8IOSCKEEN)  shown  as  solid  line.  Shaded  area 
defines  i  I  order  of  magnitude  from  the  Xu  and  Eckstein  regression  line  and 
represents  general  range  of  acceptable  values  for  dispersivity  estimates.  Note  that 
BIOSCREEN  defines  scale  as  Lp.  the  plume  Icngtii  or  distance  to  measurement  point 
in  ft,  and  employs  die  Xu  and  Eckstein  algorithm  with  a  conversion  factor  (see 
page  15). 
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DdM  Sourer;  Cielhareial.,1992 
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Figure  A. 2  Ratio  of  transverse  dispersivicy  and  vertical  dispersivity  to  longitudinal  dispersivity  data  vs.  scale 
reported  by  Gelhnr  ct  al.  (1992)  Data  includes  Golhar  s  reanalysis  of  several  dispersi/ity  studies.  Size  of  symbol 
represents  general  reliability  of  dispersivlty  estimates  Location  of  transverse  dispersivicy  r<Hationship  used  in 
BIOSCREEN  is  plotted  as  dashed  line. 
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The  Air  I'orce  vrinter  for  Environnrental  Excellence  is  distributing  BIOSCREFN  via; 


I  EPA  Center  for  Subsurface  Modeling 
I  Support  (CSMoS) 

NRMKl./SPRn 
'  r.tbBoxlldS 
I  Ada,  Cklahuma  7482M19H 


Phone:  (405)  4,36-8594 
Fax;  (405)436-8718 

Bulletin  Board;  (405)  436-8506  (14,400  baud 
8  bits  -1  stop  bit  -nu  parity). 

Web;  http;// WWW. epa, gov/ adj/l<,errlab.html 
(Electronic  manuals  will  be  in  .pdt  fo''mat;  must 
download  Adobe  A^'robat  Reader  to  lead  ami 
print  pdf  .files.) 
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APPENDIX  A.6  BIOSCREEN  EXAMPLES 


Example  1:  SWMU  66,  Keesler  AFB,  Mississippi 

•  Input  Data 

•  Fig.  1  SoLTce  Map 

•  BIOSCREEN  Modeling  Summary 

•  Fig.  2  BIOSCREEN  Input  Data 

•  Fig.  3  BIOSCREEN  Centerline  Output 

•  Fig.  4  BIOSCREEN  Array  Output 

Example  2:  LIST  Site  870,  Hill  AFB,  Utah 

•  Input  Data 

•  Fig.  5  Source  Map 

•  BIOSCREEN  Modeling  Suntmary 

•  Fig.  6  BIOSCREEN  Input  Data 

•  Fig.  7  BIOSCREEN  Ceitterline  Output 

•  Fig.  8  BIOSCREEN  Array  Output 
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BIOSCREEN  EXAMPLE  I 

Keesler  Air  Force  Base,  SWMU  66,  Mississippi 


. -y 

"ft;:  w:?; 

ikl - a.  \ ' 

Hydrugcology 

•  Hydraulic  Conductivity: 

•  Hydraulic  Gradient: 

•  Porosity: 

1.1  X  10  ’  (tin /sec) 

U.003  (ft/lt) 

0..1 

•  Slug-tests  results 

•  Static  water  level 
measurements 

•  Estimated 

Dilpersion 

Ongmal; 

•  Lon^itiidinal  Disperr.ivjiy; 

•  Traasvcrsi*  Dispersivily: 

•  Vertical  Oispersivity: 

After  Calibration: 

•  Lirngitudinai  I^ispersivity: 

•  Traasvorso  L>i  ;pcrsivity: 

•  Vertical  Dispersivity: 

13.3  (ft) 

1.3  (ft) 

0  (ft) 

32.5  (ft) 

3.25  (ft) 

0  (ft) 

•  Based  on  estimatLHl  plume 
length  of  280  ft  and 
Xu/Eckstein  relationship 

•  Based  on  calibration  to 
plume  length  (Note  this  is 
well  \  jthin  the  observed 
range  for  long,  disfx'rsivity; 
see  F'ig.  A.l  in  Appendix 

A. .3.  Remembf'r  to  convert 
from  feel  to  ini  ters  before 
using  the  chart). 

Adsorption 

•  Retardation  Fai  lor: 

•  Soil  Bulk  Density  pb: 

•  ioc: 

•  Koc: 

1.0 

1.7(kK/L) 

0.00,577„ 

B:  3»  T:  1.3,5 

E:  05  X;  240 

•  Calculair'd  from 

K  -  1  r  K(V  X  foe  X  pi  1  /n 

•  Estimated 

•  Lab  analvscs 

•  Literatii  !•  -  u.so  Koc  -  .3K 

Biodegrad  iiiion 

hlectron  Acceptor: 

Background  Cone.  fmg/L): 
Minimum  Cone  (m^/L)- 
('hange  m  Cone,  (mg/l.): 

f  lectnm  Acceptor: 

Max.  C  om .  (m)i[/L): 

Av^.  C'om .  (mg/1,): 

122  NC33  504 

2.05  0.7  2h.2 

-  0.4  •  0  ■  .t.H 

•  Based  on  March  1995 
groundwater  sampling 
program  ct>nducU“d  by 
Gnnmdwatcr  Sr’rvici's,  Inc. 

1  1.05 1  1  0.7  1  1  22.4 

lit  Uil 

36.1  7.4 

riSTl  fTT] 

Note:  BtJxinJ  valm*s  an* 

BIC>SCKLFN  input  values. 

Ge.ier.il 

•  Modeled  Are.i  Length: 

•  Modelekl  Area  Width: 

■  Simulation  Time: 

370  (ft) 

200  (ft) 

6(yr..) 

•  Basi'd  on  area  of  affected 
groundwater  plume 

•  Steady-st.ite  flow 

Source  Data 

•  Sourci'  Thickness: 

•  Skiurce  Concentration: 

in  (ft) 

(See  Figure  1) 

•  Based  on  gixdtigic  logs  and 
lumptxl  BIT'X  m'  lutormg 
tlxita 

Actual  Data 

OUTPUT 

Oi.stance  From  Source  (ft): 

BTEX  Cone.  (mg/L): 

C!enterline  (.’oncentralion: 

a)  a)  Mj  2SSJ 

5.0  1.0  0.5  O.IKIl 

•  Based  on  obseiveil 
concentrations  at  site 

StH*  F'igure  3 

Array  Conentration: 

See  Figure  4 
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BIOSCREEN  Modeling  Summary,  Keesler  Air  Force  Base,  SWMU  66,  Mississippi: 


•  BIOSCREEN  was  used  to  try  to  reproduce  tht?  movement  of  the  plume  from  1989  (the  best 
guess  for  when  the  release  occurred)  to  1995. 

•  The  soluble  mass  in  soil  and  NAPL  was  estimated  by  integrating  BTEX  soil  concentrations 
contours  mapped  as  part  of  the  site  soil  delineation  program.  An  estimated  2000  Kg  of  BTEX 
was  estimated  to  be  present  at  the  site.  This  value  represented  a  source  half-life  ot  60  years 
with  the  instantaneous  reaction  model  (the  first  value  shown  in  the  source  half-life  box  in 
Figure  2),  a  relatively  long  half-life,  so  the  2000  Kg  measured  in  1995  was  as.sumed  to  be 
representative  of  1989  conditions. 

•  The  instantaneous  reaction  model  was  used  as  the  primary  model  to  try  to  reproduce  the 
plume  length  (~  280  ft), 

•  Because  a  decaying  source  was  u.sed,  the  source  concentration  on  the  input  screen  (representing 
concentrations  6  yrs  ago)  wore  adju.stcd  so  the  .source  concentration  on  the  centerline  output 
screen  (repre.senting  concentrations  now)  were  equal  to  12  mg/L  Because  the  source  decay 
term  is  different  for  the  first  order  decay  and  instantaneous  reaction  models,  this  simulation 
focused  on  matching  the  instantaneous  reaction  model.  Tlie  final  result  was  a  source 
concentration  of  13  (i8  mg/E  in  the  center  of  the  source  7nne  (note  on  tlie  centerline  output  the 
source  concentration  Is  12.021  mg/L). 

•  The  initial  run  of  the  instantaneous  reaction  model  indicated  that  the  plume  was  too  long. 
Thi.s  indicates  that  there  is  more  mixing  of  hydrocarbon  and  electron  acceptors  at  the  site 
than  is  predicted  by  the  model.  Therefore  the  longitudinal  dispersivity  was  adjusted 
upwards  (more  mixing)  until  BIOSCREEN  matched  the  observed  plume  length.  The  final 
longitudinal  dispersivity  was  32.5  ft. 

•  Asa  check  the  first-order  decay  model  was  used  witli  the  BIOSCREEN  default  value  of  2 
yrs.  This  run  greatly  overestimated  the  plume  length,  so  the  amourt  of  biodegradation  was 
increased  by  decreasing  the  solute  half-life.  A  good  match  of  the  plume  was  reached  with  a 
.solute  half-life  of  O.l.S  years. 

•  As  shown  in  Figure  3,  BIOSCKEEN  matches  the  observed  plume  fairly  well.  The 
instantaneous  model  is  more  accurate  near  the  source  while  the  first  order  decay  model  is 
more  accurate  near  the  middle  of  the  plume.  Both  models  reproduce  the  actual  plume 
length  relatively  well. 

•  As  shown  in  t'igiire  4,  the  current  plume  is  estim.ited  to  contain  7.8  kg  of  BTEX.  BIO.SCREHN 
indicates  that  the  plume  under  a  nivdegradation  scenario  would  contain  126.3  kg  BTFiX,  In 
other  words  BIOSCKEEN  indicates  that  94'Xi  of  the  BTEX  mas.s  that  has  left  the  source  since 
1989  ha.s  biodegraded. 

•  Most  of  the  .source  mass  postulated  to  be  in  place  in  1989  is  .still  there  in  1996  (2(X.X)  kg  vs.  1837 
kg,  or  92%  left). 

•  The  current  plume  contains  l.Oac  ft  of  contaminated  water,  v/ith  1.01'^  acre-ft/yr  of  water 
being  contaminated  as  it  flows  through  the  soums  Because  the  plume  is  almost  at  steady 
state,  1.019  ac-ft  of  water  become  contaminated  per  year  with  the  same  amount  being 
remedi.ited  every  year  due  to  in-situ  biodegradation  and  other  attenuation  processes. 
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DISSOLVED  HYDROCARBON  CONCENTRATION  ALONG  PLUME  CENTERUNE  (mftil.  *tZ-0) 


Figure  3.  Centi^rline  Output.  Keesler  Air  Force  Base,  Mississippi. 


DISSOLVED  HYDROCARBON  CONCENTRATIONS  IN  PLUME  (rag/L  atZ=0) 
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Figure  4.  Conce-i'-ration  Outp-c.  Keesier  Ai'-  Force  Base.  Miss'ss  pp 
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EXAMPLE  2 

Hiii  Air  Force  UST  Site  870,  Utah 


DATA^YFT-lg-S 

Hydrogaolocy 

•  Hydraulic  Conductivity; 

•  Hydraulic  Gradient: 

•  Porosity: 

8.05  X  10'  (cm/sec) 

0.048  (It/ll) 

0.25 

•  SluR-tests  results 

•  Static  water  level 
measurements 

•  Estimated 

Dispcnion 

Ongirul 

•  LaOngitudinal  Dispersivity: 

•  Transverse  Dispersivity: 

•  Vertical  Dispersivity: 

28.5  (ft) 

2.85  (ft) 

It  (ft) 

•  Ba*^  on  estimatkxl  plume 
length  of  1450  ft  and  Xu's 
disj>ers»vilY  h^rmula 

•  Note:  No  calibration  was 
nLxressary  to  match  the 
ob.sorveii  plume  length. 

Adsorption 

•  Retardation  Fa(.-tor: 

•  tjoil  Bulk  Dc'asitv  pb: 

•  fix': 

•  Koc: 

1.3 

1.7(kK/') 

0.08% 

1):  38  T:  135 

K:  95  X  240 

•  Cak'ulati'd  from 

K  1  ♦Koc  X  hv  X  pb/n 

•  Estimated 

•  l^h  analysis 

•  Literature  -  use  Kih*  35 

Biodegradation 

lili’ctron  Acceptor: 
liackKround  Cone. 

Minimum  (ioiic  (rnn/L.): 
Change  in  Cone.  (m>;/l.): 

lilectron  Acceptor: 

M.i\.  Cone.  (mK/1.): 

Avk-  Cone.  (m>;/l.): 

8.0  17.0  KXi 

-  0.22  ■  0  ■  0 

1  .5.78 1  1  17.1)1  1  lOit'l 

li:  cm 

•4).S  2.04 

□13  iMlil 

Note:  Bo^ed  values  are  DiOSC'KTL'N 
input  values. 

•  Based  n  lulv  1994 
(groundwater  samplin(> 
program  conducted  by 
Parsons  Enniius'rinj; 
Science,  Ine. 

General 

•  Mliuficd  Ar<*a  Length 

•  MiKJeled  Area  Width: 

•  Smuilalion  I  ime: 

1450  (tt) 

,321)  (ft) 

5  (vrs) 

•  Basr  d  on  area  of  .iffe(  led 
groundwater  f>hime 

•  Steady-stale  flow 

Source  Data 

•  Sourer  Thickness' 

•  Source  Concentration: 

10  (ft) 

{StH*  Figure  S) 

•  Basis.!  on  gisilogu’  l(»gs  and 
liini[)i’d  BTLX  monitoring 
data 

Actual  Data 

nistanee  Irom  Source  (It): 

1)  I  PX  Cone.  (niK/I.): 

m  im  1350  mj 

K.O  III  11.1)2  O.tKtfi 

*•  IJiiM*d  on  obsi*rvevI 

concentration  contour  .it 
site  (sis'  Figure  .5) 

OUTPUT 

Centerline  ConcGntraiior.. 

rigim*  ~f 

Array  (  onceritration: 

Sev  ITgiiri*  H 

- —  — 1 
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Source^Zone^ssumptiof^ 


um-i 
K4WS  f.m 


Source  of  Data:  Wiedemeler  et  al ,  1P95b 


LEGEND 


'  Monitoring  vyell  location 
■  July  1994  Geoprotie  sampling  location 
SO-  BTEX  concentration  Isopleth,  mg/L,  July  1994 
Affected  Soil  Zone 


BIOSCREEN  SOURl  F  ZONE 
ASSUMPTIONS 

U5T  Sile  070,  Hill  AF-B,  Utah 

FIGURE  5 
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BIOSCREEN  Modeling  Summary  Hill  Air  Force  Base,  UST  Site  870,  Utah: 


•  BIOSCREEN  was  used  to  try  to  reproduce  the  movement  of  the  plume. 

•  An  infinite  source  was  assumed  to  simplify  the  modeling  scenario.  The  source  was  a.ssumed  to 
be-  in  the  high  concentration  zone  of  the  plume  area  (see  Figure  S).  Note  that  the  zom“  of 
affected  soil  was  quite  large;  however  much  of  the  affected  soil  zone  downgradient  of  tlu“ 
.source  wa.s  relatively  low  concentration. 

Two  modeling  approaches  could  be-  applied:  1)  a.ssuming  the  .source  zone  is  just  downgradient 
of  the  affected  soil  area  (near  well  EI’A-82-C)  and  ignoring  the  area  upgradient  ot  the  this 
point,  ind  2)  modelmg  most  of  the  plume  with  souav  near  MW-1.  Alternative  1  is 
theoretically  morx’  accurate,  as  HIOSCKEEN  cannot  aevount  for  the  contributions  tnim  any 
affected  soil  zone  downgradient  of  the  .source.  At  the  case  of  Hill  AFB,  however,  it  was 
a.ssumed  that  the  contributions  from  this  downgradient  affected  soil  were  relatively  muior 
and  that  the  main  priK'ess  of  interest  was  the  length  of  the  plume  from  the  high- 
concentration  source  zone.  Therefore  Alternative  2  was  modeled,  with  the  note  that  the 
middle  ot  the  actual  plume  may  actually  have  higher  concentrations  than  would  Ix’ 
expected  due  to  th  contaminants  in  the  downgradient  affected  soil  z,on*‘. 

•  The  instantaneous  reaction  model  was  ixsed  as  the  primary  mixiel  to  try  to  nqinxluce  the 
plume  length  (-  280  ft)  as  shown  in  I'lgure  7. 

•  The  initial  run  ot  the  instantaneous  reaction  mrxJel  repnxJucixi  the  existing  plume  without 
.my  need  for  calibration  of  dispersivity. 

•  As  a  check  the  first-orde,'  decay  nuxlel  w.is  usixi  with  the  BIOSCREEN  default  value  of  2 
yrs.  This  run  greatly  overestimated  the  plume  length,  .so  the  amount  of  biodegradation  w.is 
increa.sed  by  decreasing  the  .solute  hall-life.  A  h.ilf-life  value  of  0.1  years  w.is  required  to 
match  the  |)lume  length,  .ilthough  the  m.itch  in  the  middle  in  the  plume  was  much  poorer. 

•  As  shown  in  Figure  7,  BIOSC’REEN  matches  the  ob.st'rved  plume  f.iirly  well  I'he 
inst.intaiu'ous  nuxlel  is  more  accurate  near  the  soutvv  while  the  first  order  dec.iy  model  is 
inoR'  accurate  near  the  iiiiddk'  of  the  plume.  Both  nuxlels  lepnxliu’e  the  .letii.il  plume 
li’iigth  rel.ifively  well. 

As  shown  in  Figure  H,  the  model  w.is  iituible  to  ealcul.ite  the  ni.i.ss  h.il.inces.  A  quick 
i-valuatuui  shows  the  reason:  witli  .i  ser*p.ige  velix  ity  ot  IN)*'  ft/yr  and  .i  5  ye.ir  sinuil.ition 
time,  the  iindegr.uled  plume  should  hi-  over  8<KK)  ft  long.  Iks  .luse  the  ni.l.ss  b.il.iiice  is  h.ised 
on  a  coiiijuinson  ol  ,i  complete  undegraded  plume  vs.  a  degr.uled  pliimi-,  .i  modi'l  .ire.i  leiigtli 
ot  8(K)0  tl  would  Ix'  requirr'd  tor  HiDSCUEEN  to  compk-te  the  iiia.ss  li.il.ince  ealeiilatum. 
I  heretore  two  nins  Would  Ix' lUH'ded  to  complete  the  simulation:  1)  ,i  nni  with  .i  modeled 
length  ot  I'lSO  ler-l  to  c. ilibr.il.'  .md  <'v.tlute  fbe  mateli  to  existing  data,  and  2)  .i  nin  with  .i 
modeled  length  of  SiXMl  tt  to  do  the  ni.l.ss  b.ilaiiee.  The  results  ol  tlie  srsoiul  nm  (ch.inge  ot 
modi’l  area  length  from  14.‘i0  ft  to  HOfK)  ft)  indicate  tliat  over  W'/,  of  the  m.i.ss  tli.il  has  left 
the  source  has  biodegraded  by  tlie  time  gronndw.iter  has  traveled  HSO  ft. 
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